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ABSTRACT 


The  interaction  between  a  grid  modulated  electron  beam  of  finite 
diameter  and  a  stationary  plasma  col\imn,  in  the  absence  of  a  magnetic 
field,  was  investigated  both  theoretically  and  experimentally. 

Experiments  performed  with  a  beam-plasma  tube,  with  modulation  of 

the  electron  beam  by  grids,  revealed  a  maximum  interaction  when  the 

modulation  frequency  of  the  beam  was  equal  to  or  near  the  plasma  frequency 

of  the  plasma.  An  electronic  gain  of  7*5  db  due  to  the  plasma  was 

observed.  The  output  povrers  v^ere  low  with  no  net  gain  in  any  condition, 

due  to  the  low  impedances  of  the  cavities  which  were  used,  and  to  the 

low  beam  current  available.  Nevertheless,  it  v/as  possible,  by  driving 

the  beam  to  saturation  power  levels,  to  obtain  fundamental  components  of 

the  beam  current  such  that  *  ^*5  near  the  resonance  condition. 

With  only  the  electron  beam  on  the  plasma  being  off,  the  saturation 

output  current  measured  was  about  •  These  large  saturation 

currents  which  were  observed  may  be  attributed  to  the  fact  that  when 

a>  <  CD  the  plasma  behaves  like  an  inductive  medium  which  has  a  similar 
PP 

effect  to  the  inductively  detuned  intermediate  cavity  in  a  three-cavity 
klystron. 

Second  harmonic  signals  \;ere  detected,  and  a  maximum  peak  of  the 

output  signal  was  observed  v/hen  the  plasma  frequency  of  the  plasma  was 

twice  the  modulation  frequency  of  the  beam.  Another  smaller  peak  was 

observed  at  f  ~  f  .  The  second-harmonic  currents  were  shown  to  increase 
P 

approximately  by  a  factor  of  two,  at  input  saturation  powers,  when  the 
plasma  density  was  increased  frcxn  zero  to  a  value  such  that  f ^  2f  . 

The  theoretical  model  which  was  the  subject  of  analysis  consisted  of  a 
perfectly  conducting  cylindrical  waveguide  completely  filled  by  a  loss¬ 
less,  uniform  plasma.  The  electron  beam  was  assumed  to  be  modulated  by 
a  gridded  gap. 

Two  pairs  of  field  solutions  inside  the  beam  region  are  shown  to 
be  possible.  One  pair  of  solutions,  not  excited  with  a  gridless  gap 
modulation,  was  associated  with  space-charge  density  variations  inside 
the  beam.  The  propagation  constants  of  these  volume  space-charge  waves 
v;ere  the  seune  as  in  the  infinite  beam  plasma  system.  The  other  pair  of 
solutions  for  the  electric  fields  within  the  beam  corresponds  to  the 
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condition  for  which  there  is  no  charge  accumulation  in  the  beam.  These 
solutions  are  associated  with  charge  rippling  at  the  beam-plasma  interface. 
The  propagation  constants  of  these  surface  waves  depend  on  the  transverse 
dimensions  of  the  beam  and  plasma.  These  svirface  waves  grow  exponentially 
with  distance,  when  o)  <  ,  at  a  rate  which  is  smaller  than  the  growth 

rate  of  the  volume  space-charge  waves. 

It  was  shown  also  that,  under  conditions  for  existence  of  growing 
waves  in  the  system,  strong  harmonic  generation  becomes  possible  as  a 
result  of  modifications  of  space-charge  forces  on  the  beam,  due  to  the 
presence  of  the  plasma. 

Expressions  for  the  second-harmonic  conppnents  of  the  relevant  rf 
quantities  in  the  infinite  beam-plasma  system  were  derived,  indicating 
the  existence  of  growing  second-harmonics. 
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CHAPTER  I 
INTRODUCTION 

The  purpose  of  this  study  was  to  investigate  the  interaction  effects 
between  a  velocity-modulated  electron  beam  and  plasma,  effects  which  result 
in  amplification  and  harmonic  generation.  Tae  existence  of  spatially 
growing  waves  was  predicted  some  years  ago  by  Bohm  and  Gross. ^  Experiments 

2  3 

performed  in  1958  by  Boyd,  Gould  and  Field,  and  in  Russia  by  Bogdanov*^ 

et.  al.,  successfully  demonstrated  the  existence  of  such  growing  waves. 

More  recently,  similar  experiments  using  a  thermally  generated  cesium 

4 

plasma  were  reported  by  Allen  and  Kino.  In  these  experiments,  a  velocity- 
modulated  electron  beam  was  allowed  to  transverse  a  plasma  region,  and 
maximum  interaction  between  the  beam  and  the  plasma  was  obtained  at  a 
modulation  frequency  which  was  equal  to  or  near  to  the  plasma  frequency. 

We  shall  be  primarily  concerned  with  the  amplification  and  the  associated 
harmonic  generation  mechanisms,  which  we  snow  to  be  related  to  tne  velocity 
modulation  process  and  subsequent  buncning  of  electrons  of  t  le  beam,  in  the 
plasma  medium. 

Chapter  II  is  devoted  to  an  analysis  which  makes  use  of  the  Lagrangian 
approach  to  determine  the  beam-plasma  interaction  mechanism.  It  is  shown 
that  a  high  harmonic  content  would  be  expected  to  exist  in  the  current  in 
the  bunched  beam  supporting  a  growing  space-charge  wave.  Only  the  case  of 
a  plasma  and  electron  beam  of  infinite  extent  is  considered  in  this  chapter. 

Chapter  III  deals  with  a  more  practical  configuration  which  consists 
of  an  electron  beam  of  finite  diameter  passing  along  the  axis  of  a  plasma 
column,  with  no  external  dc  nagnetic  field  applied  and  velocity-modulated 
by  grids. 

It  is  shown  in  the  analysis  of  this  problem  that,  with  modulation  by 
grids,  the  space-charge  waves  proi>agating  along  the  beam  increase  ex¬ 
ponentially  with  distance  at  a  rate  of  growth  which  is  higher  than  the 
one  corresponding  to  the  wave  solution  in  the  case  of  gap  or  helix  mod¬ 
ulation.  The  excitation  of  the  individual  space-charge  waves  by  a  gridded 
gap  is  determined  and  an  expression  for  the  rf  beam  current  at  the  out¬ 
put  grid  is  derived.  These  solutions  which  lead  to  these  high  growth 
rates  correspond  to  the  case  in  which  there  is  a  charge  bunching  within 

the  beam;  they  are  nonsolenoidal  and  p  /  0  inside  this  region. 
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When  the  modulation  frequency  of  the  electron  beam  is  above  the  plasma 
frequency  of  the  plasma,  it  is  shown  in  the  text  that  there  is  no  field  in 
the  plasma  outside  the  beam  corresponding  to  these  nonsolenoidal  waves. 

When  the  o  <  the  propagation  constants  become  complex  and 

the  proof  given  in  the  text  does  not  apply#  It  can  be  shown  however  by 
an  analagous  treatment  of  Bessel  functions  of  complex  argument  that  the 
fields  of  the  nonsolenoidal  waves  vanish  in  this  case  also,  ^ 

The  method  which  is  used  throughout  this  analysis  is  the  superposition 
method,  in  contrast  to  the  Bblti^niann  approach  or  the  Fourier  analysis  of 

7 

density  distrioution.  The  reader  is  referred  to  Buneman  for  a  discussion 
of  the  latter  two  methods.  According  to  the  superposition  approach  we 
treat  the  streams  of  electrons  (beam  and  plasma)  by  a  perturbation  method 
leadin.,  to  the  Eulerian,  small-signal  description  of  motion  of  electrons. 

The  term  plasma  is  used  in  this  report  to  denote  an  assembly  of  equal 
numbers  of  electrons  and  positive  ions  and  perhaps  neutral  molecules. 

The  ions  are  assumed  stationary.  Unless  otherwise  specified,  thermal 
motions  and  collisions  are  neglected. 

The  nonlinear  aspects  of  the  beam-plasma  interaction  problem  are 
investigated  in  Chapter  IV  by  means  of  a  complete  formal  deri\ration  of 
the  relevant  second  harmonic  quantities. 

The  experimental*  work  is  described  in  the  final  chapter.  The 
purpose  of  the  experiment  was  to  investigate  the  beam-plasma  interaction 
near  the  plasma  frequency  of  the  plasma,  under  the  conditions  of  tae 
analysis  made  in  Chapter  III,  i.e.,  (a)  no  magnetic  field,  and,  (b)  mod¬ 
ulation  by  grids.  Other  interesting  effects  concerned  with  beam  saturation 
data  and  second  harmonic  generation  are  also  reported  and  discussed. 
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CHAPTER  II 


ELECTRON  BEAM-PLAEMA  INTERACTION  MECHANISM;  LAGRANGIAN  APPROACH;  HARMONIC 

GENERATION 


2.1  INTRODUCTION 

Before  undertaking  (in  the  next  chapter)  the  analysis  of  the  inter¬ 
action  of  a  finite  electron  beam  with  a  plasma  column  through  which  it  passes, 
it  is  of  interest  to  coAsider  the  essential  aspects  of  the  beam  plasma 
interaction  which  are  apparent  in  the  simple  one-dimensional  treatment 
of  the  problem. 

In  this  introductory  chapter  we  shall  deal  then  with  the  ideal 
situation  of  a  one-dimensional,  infinite  cross-section  beam  model  and 
an  unbounded  cold  plasma. 

Our  analysis  shall  be  restricted  to  small  signals.  The  possibility 
of  cross-over  of  beam  electrons  is  excluded,  so  that  a  unique  velocity 
vector  is  associated  with  any  point  of  the  beam  at  a  given  time.  It  is 
shown  that  each  infinite  sheet  of  electrons  constituting  the  plasma  will 
oscillate  about  its  equilibrium  position,  provided  that  the  sheet  does 
not  cross  another  sheet  in  the  course  of  its  motion. 

It  is  the  purpose  of  this  chapter  to  investigate,  under  these 
circumstances  and  without  going  into  details,  the  nature  of  the  space- 
charge  effects  of  a  modulated  electron  beam  drifting  along  the  z-dlrection 
through  the  plasma  medium.  To  study  the  influence  of  the  plasma  upon 
the  space-charge  waves  on  the  moving  beam,  we  chose  the  Lagrangian 
description  where  attention  is  primarily  concentrated  on  the  trajectories 
of  individual  electrons.  It  is  indicated,  in  this  elementary  analysis, 
that  the  plasma  affects  the  bunching  in  the  beam  in  such  a  way  that, 
under  some  conditions,  growing  space-charge  waves  are  supported  by  the 
electron  beam.  Also,  under  these  same  conditions,  strong  harmonic 
generation  might  be  possible  as  a  result  of  modifications  of  space-charge 
forces  on  the  beam,  due  to  the  presence  of  the  plasma. 

2.2  EQUATIONS  FOR  THE  UNCOUPLED  ELECTRON  BEAM  AND  PLASMA 

The  system  of  variables  useful  in  the  small  signal  analysis  and 

8 

suitable  for  our  present  purposes  is  described  by  Bobroff  es  the 
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"polarization  system"  which  In  the  absence  of  a  dc  velocity  (as  In  our 
case  of  a  stationary  plasma)  becomes  Identical  to  the  Lagranglsji  system. 

In  this  section  a  brief  discussion  of  the  dependent  variables  which 
play  Important  roles  In  the  analysis  which  follows  Is  presented.  Reference 
should  be  made  to  the  Bobroff  paper  for  further  details. 

We  shall  first  define  the  ac  position  vector  or  displacement  vector. 

It  Is  known  that  when  modulation  Is  applied  to  the  Input  gap  the  electrons 
constituting  the  beam  are  displaced  from  their  dc  positions.  This  dis¬ 
placement  will  be  called  the  displacement  vector  and  will  be  denoted  by 
.  To  first  order  It  Is  a  function  of  the  dc  position  and  also  an 
explicit  function  of  the  time.  Thus  we  write 

^^^Ob  "  “  *0b  ^b^^OL  ^ 

where  z^  denotes  the  dc  position  of  the  beam  electrons. 

The  "polarization"  or  electric  displacement  arises  from  the  charge 
separation  caused  by  the  displacement  of  electrons  from  their  equilibrium 
positions.  For  the  electron  beam.  It  Is  defined  to  first  order  In 
terms  of  the  above  displacement  as 

\  -  Pob  •  tb  (2-2) 

where  Is  the  dc  charge  penalty  of  the  beam.  The  minus  sign  In  this 

expression  accounts  for  the  fact  that  we  deal  with  displacements  of 
negatively  charged  sheets. 

The  other  first-order  dependent  variable  Is  the  ac  velocity  which 
Is  denoted  by  v.  . 

D 

The  Eulerlan  charge  and  linearized  current  density  may  be  expressed, 
as  Is  shown  by  Bobroff,  in  terms  of  the  above  variables.  Ihe  ac  cheurge 
density  of  the  electron  beam  is 


The  linearised  current  density  of  the  becm  is  given  by  the  relation 

\  ^  •  ^b^  Vx  (pQb^b 

It  can  be  shown  (see  reference  C)  that  the  follov/ing  relation  holds 

between  v,  and  : 

°i  ^ 


(2.5) 


Similarly,  the  equation  of  motion  for  the  beam  electrons  in  this 
system  of  variables  is  to  first-order: 


O)  , 

=  -EL. 


Dv. 


Ob 


(2.0) 


where 


O) 


pb 


'Z  ”‘o 


(2.7) 


is  the  pli’sma  frequency  of  the  electron  beam. 

It  follov.'G  from  Eqs.  (2.2),  (2-5)  and  (2.6)  that 


3  c)  2 

ySt  *  “0  •  "  ; 


2 

5  =  -2L.  D 

^b  •  ^b 

"ob 


(2.8) 


In  a  frame  of  reference  moving  at  the  dc  beam  velocity,  the  luantities 
5^  and  obey  the  following  differential  equations: 


(2.9) 


where  t  * 


represents  time  in  the  moving  system. 


The  solution  of  Eq.  (2.8),  which  gives  the  displacements  of 

the  beam  electrons  drifting  along  the  z  direction  with  a  dc  velocity 

9 

Uq  ,  Is  of  the  form 

sln(ajt  -  slnp^^z^  (2.10) 


where 


A 


V 


20) 


pb 


with 


rf  modulation  voltage 
Vq  dc  beam  voltage 


and  the  denotes  the  dc  equilibrium  position  of  the  beam  electrons, 
and  it  is  assumed  that  at  t  =  t^  we  have 

=  ° 


z  =  Uq(1  +  I  slnojt^)  ;  a  «  1  . 

It  follows  similarly  that  the  equation  of  motion  for  the  plasma  electrons 
is 


5t 


(2.11) 


where 


(U 


^  /!2e! 

V“^o 


designates  the  plasma  frequency  of  the  plasma,  and 

V 

^1  dr 

is  the  ac  velocity  of  the  plasma  electrons. 


(2.12) 
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By  coinblning  together  Eqs.  (2.11)  and  (2.12),  and  by  using  the 

relation  D  =  -  p-  £  ,  we  obtain  for  the  plasma  small  amplitude  equations 

P  Op  P 

of  simple  harmonic  motion: 


+  05 


PP 


0 


+  05 


PP 


D 

P 


=  0 


(2.13) 


These  equations  represent  plasma  oscillations  which,  In  the  case  we 
are  considering  (a  stationary  plasma),  do  not  propagate,  since  any 
disturbance  would  in  the  absence  of  collisions  or  thermal  velocities 
persist  indefinitely  at  the  location  of  the  disturbance. 

2.3  the  inhcmogeneous  differential  equations  for  the  linear  interaction 

OF  A  BEAM-PLASMA  SYSTEM' 

We  have  been  considering  In  the  previous  section  the  electron  beam 
and  the  plasma  as  two  Independent  uncoupled  systems.  We  shall  be  con¬ 
cerned  now  with  the  Interaction  that  tedces  place  between  the  plasma  and 
the  beam  when  they  are  coupled  together. 

We  shall  restrict  ovirselves  to  linear  Interaction.  This  means  that 
we  shall  use  only  the  linearized  equations  of  motion  which  were  derived 
In  section  2.2. 

nie  displacements  of  the  electrons  of  the  beam  and  plasma  axe  changed 
by  the  coupling.  These  displacements  are  therefore  denoted  by  and 

5^  ,  respectively. 

The  electric  field  In  the  beam-plasma  system  is 


where 


The  equation  of  motion  for  the  plasma  electrons  is 


(2.15) 


For  a  one-dlmenslonal  system,  the  total  longitudinal  current  (con¬ 
vection  plus  displacement  current)  must  he  zero.  It  follows  that 


(2.16) 


Substitution  of  Eq.  (2.^)  in  Eq.  (2.l6)  yields  the  result 


d  d 

dt  ^  °  Jt  Op  p 


(2.17) 


since  the  second  term  in  Eq.  (2.4)  is  zero  for  4he  one  dimensional  hesB. 

Partial  differentiation  with  respect  to  time  of  the  above  equation 
gives  the  following  expression: 


b^i. 

i?* "  ■  ■  j?"  ■  "Of  ■  j? 


(2.18) 


into  which  may  be  substituted  Eq.  (2.15)  to  yield  the  inhonogeneous 
differential  equation  for  : 


+  (D 


PP 


(2.19) 


The  equation  of  motion  for  the  beam  electrons  is 


(2.20) 
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This  equation,  together  with  the  preceding  one,  constitutes  a  system 
of  equations  which  may  then  he  solved  by  assuming  that  all  quantities  vary 
as  exp  (cut  -  pz).  The  following  dispersion  relation  is  then  obtained 
from  Eqs.  (2.19)  and  (2.20): 


CO  2 

-  -22_  . 

2 

CO 


CO 


£L 


(co  - 


(2.21) 


from  which  it  may  be  seen  that 


P  =  Pg  ±  h  , 


h  = 


(2.22) 


where 


CO  CO  . 

u^ 

0  0 


This  equation  indicates  that,  when  the  electron  beam  is  modulated  at  a 

frequency  co  above  the  plasma  frequency  of  the  plasma  co^^  ,  there  are 

two  propagating  waves  that  can  be  excited:  one  with  a  phase  velocity 

greater  than  u^  ,  and  the  other  with  a  phase  velocity  less  than  u^  . 

When  CO  <  CO  ,  the  denominator  of  the  right  hand  side  of  £q.  (2.22)  is 
PP 

negative,  ^  becomes  complex  and  two  waves  are  found:  one  growing  with 
distance,  and  the  other  decreasing  with  distance.  By  introducing  Eq. 

O 

(2.14)  into  Eq.  (2.19),  and  writing  Pq^/Pq^  =  (“plj/^pp)  >  following 
differential  equation  for  the  displacements  of  the  electrons  constituting 
the  plasma  is  obtained: 


+  CO 


PP 


(2.23) 


-  9  - 


The  right-hand  side  of  this  equation  can  be  written  in  terms  of  the  beam 
current  density,  yielding 


+  05 


PP 


-  ^  . 


cue. 


Instead  of  free^  undamped  oscillations  of  the  plasma  electrons  as 
given  by  Eq.  (2.13)  in  the  presence  of  the  electron  bean^  we  now  have 
forced  oscillations  in  the  plasma  with  the  driving  force  being  due  to 
the  ac  modulated  beam. 

By  assuming  as  before  that  the  displacements  of  the  plasma  electrons 
varies  as  e*^^  ,  it  will  be  seen  that 


0) 


'p* 


£L 


O)  «  0) 


PP 


(2.24) 


This  equation  established  the  relation  between  the  steady-state 
displacements  of  the  electrons  which  constitute,  respectively,  the 
electron  beam  and  the  plasma. 

From  Eqs.  (2.14)  and  (2.24)  it  is  seen  that  the  electric  field 
originating  from  this  displacement  of  charges  is 


“l  ■  -  '’Ob 


^b  "  ^Op^p 


-  P 


Ob 


05 

0)  > 


•  ^b 


(2.25) 


Ilie  resulting  force  on  a  beeun  electron  due  to  this  displacement 
field  is 


-  ep 


F  =  —  D, 


•  «b 


(2.26) 


On  examination  of  this  equation,  it  may  be  seen  that  (F/|^) 
reverses  sign  as  the  frequency  is  changed  from  05  >  05  to  05  <  05 


PP 


PP 
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Consequently  the  displacement  field  bunches  the  beam  when  oo  <  , 

Instead  of  debunching  it  as  it  wo\;ild  be  if  there  were  no  plasma  present 
or  if  0)  >  CDpp  .  We  should  expect  that  because  of  this  tight  bunching 
of  the  beam  electrons,  there  would  be  a  high  harmonic  content  in  the 
resulting  current. 

Equation  (2.26)  can  also  be  written  in  the  following  form: 


F 


with  the  result  that 


^2 


.2 


(2.27) 


(2.28) 


in  a  coordinate  system  moving  with  the  beam. 

(knparing  the  above  expression  with  Eq.  (2.9)>vhich  was  derived  in 
the  absence  of  plasma,  gives  a  simple  interpretation  to  the  effects  that 
the  plasma  has  upon  the  motion  of  the  beam  electrons.  We  Interpret  the 
presence  of  plasma  as  having  the  effect  of  changing  the  plasma  frequency 
of  the  electron  beam  from  its  value  ,  in  the  absence  of  plasma,  to 

the  effective  value 


0) 


qb 


(2.29) 


In  a  system  of  reference  stationary  with  respect  to  the  beam,  Eq. 
(2.28)  becomes 


(2.30) 


11 


After  expanding  the  left  hand  side  of  the  above  differential  equation, 
we  find  that 


bz^  Uq  3t5z  St^ 


(2.31) 


This  second-order  homogeneous  differential  eqviatlon  will  be  derived  again, 
in  Chapter  IV,  by  using  the  £ulerian  approach,  as  an  introduction  to  the 
nonlinear  analysis  which  will  be  developed  in  that  Chapter. 

2.k  ELECTRON  CROSSOVER;  NONLINEARITY 


The  fact  that  the  amplitude  of  ,  in  Eq.  (2.2U)  for  example, 
becomes  infinite  for  a>  =  illustrates  the  phenomenon  of  resonance. 

Care  should  be  taken,  however,  in  discussingthe  implications  of  this 
condition,  since  our  linear  differential  equations  hold  only  for  sinsdl 
oscillations  of  electrons  within  the  framework  of  the  small-signal  theory 
of  space.<charge  debunchlng  forces. 

As  pointed  out  before,  when  the  modulating  frequency  f  is  below 
fp  and  approaches  f^  ,  there  exist  strong  space-charge  forces  on  the 
beam  electrons  in  such  a  direction  to  reinforce  the  bunching  process  set 
up  in  the  beam  by  velocity  modulation. 

Under  this  condition,  the  phenomenon  of  electron  crossover  would  be 
expected  to  occur.  It  is  known  fron  the  study  of  the  bvinchlng  motion  in 
klystron  theory  that  crossover  of  electrons  exists  whenever  the  following 
inequ8^.ity  holds 


1  + 


(2.32) 


The  displacements  of  the  plasma  electrons  were  shown  to  be  related  to 
the  beam  electron  displacements  according  to  the  equation 
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It  follows  then  that 


By  assuming,  as  is  usual,  that 


_h_ 


< 


1 


9 


(2.33) 


(2.34) 


we  have,  from  Eq.  (2.33)  and  the  crossover  criterion  [Eq.  (2.32)]  the 
following  result: 


^6' 

^6’ 

b 

> 

_E 

bz 

dz 

indicating  that  the  beam  electrons  cross  over  before  such  a  phenomenon 
occurs  for  the  plasma  electrons. 

The  relationship  between  the  ac  charge  densities  of  the  electron 
beam  and  plasma  may  be  easily  derived.  We  know  that 


Pb 


1 


(a.35a) 


The  following  equation  is  then  obtained,  with  the  help  of  Eq.  (S.33): 


(U 


2 
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This  equation  can  be  written  in  the  following  form: 


By  using  some  typical  experimental  figures  it  may  be  seen  that 


In  our  experiment  it  will  be  shown  that  this  qmntity  is  of  the  order 

+4 

of  10  .  Therefore,  we  may  have 


Furthermore,  the  same  equation  (2.2^)  may  be  used  to  derive  a  relation 
between  the  plasma  velocity  and  the  ac  beam  velocity.  We  obtain 

We  know  however  that 


It  results  in  the  following  expression  for  the  plasma  velocity  in  terms 
of  the  ac  charge  density  and  ac  velocity  of  the  electron  beam: 


It  is  apparent  then  that  we  may  have 
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These  inequalities  for  relative  magnitudes  of  charge  densities  and 

velocities  do  not  determine  the  relative  meignitudes  of  p  ,v  .  and 

P)1  PA 


Of 


However  a  more  careful  consideration  of  common  values 

.2 


and 


(a>  /a>)‘ 

pp/  ' 


woxild  indicate  that  in  the  Interesting  range 


of  behavior  one  would  usually  have 


P  V 

Pi  Pi 


Pb7' 


1  ^1 


This  means  that  second-order  nonlinear  effects  may  become  first  more 
important  in  the  plasma.  In  a  nonlinear  analysis  to  be  carried  out  in 
Chapter  IV,  the  nonlinear  term  (p,  ^v  ,  )  which  enters  in  the  definition 
of  the  second-order  plasma  current  density  may  be  considered  a  more 
relevant  term  than  the  corresponding  second-order  term  (p^  ^v^  in 
the  expression  for  the  beam  current  density. 
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CHAPTER  III 


INTERACTION  OF  A  GRID-MODULATED  ELECTRON  BEAM  WITH  A  PLASMA  (NO  EXTERNAL 

MAGNETIC  FIELD) 


3.1  INTRODUCTION 

This  chapter  is  concerned  with  an  .analysis  of  the  interaction  of 
an  electron  beam  of  finite  radius  with  a  plasma.  More  specifically,  as 
shown  in  Fig.  1,  we  consider  a  perfectly  conducting  cylindrical  waveguide 
of  diameter  2b  completely  filled  with  an  ideal  electron  plasma.  An 
electron  beam  of  diameter  2a(a  <  b)  drifts  with  a  constant  uniform 
velocity  Uq(Uq  «  c)  through  this  plasma  column  after  being  velocity- 
modulated  by  ideal  grids.  Identical  output  grids  are  used  as  parts  of  an 
external  cavity  where  the  demodulation  process  takes  place.  It  is  assvmed 
that  there  is  no  dc  magnetic  field  applied  and  the  electron  beam  is 
focused  by  positive  ions. 

The  problem  is  to  obtain  solutions  for  the  space-charge  waves  in 
this  system,  to  combine  these  solutions  to  match  the  boundary  conditions 
at  the  input,  and  to  derive  expressions  for  the  rf  beam  current  at  the 
output  grids.  Feenberg^  developed  the  theory  of  small  signal  bunching  In 
a  beam  of  finite  cross-section  In  a  drift  tube,  with  positive  ion  focusing 
and  zero  dc  magnetic  field.  The  following  analysis,  where  modulation  is 
assumed  to  be  by  grids,  is  based  on  Feenberg's  work  which  is  modified  to 
take  into  account  the  presence  of  plasma  in  the  drift  space. 

It  will  be  shown  that  the  plasma  has  the  effect  of  changing  the 
expressions  for  the  field  quantities  in  such  a  way  that,  under  certain 
conditions,  waves  increasing  exponentially  with  distance  become  possible. 
Moreover,  with  modulation  by  grids,  volume  space-charge  waves  involving 
charge  density  variations  are  shown  to  be  possible  solutions  of  the 
problem.  It  is  important  to  note  that  their  growth  rate  is  greater  than 
the  corresponding  growth  rate  for  the  wave  solutions  in  the  case  of 
modulation  of  the  electron  beam  by  a  gap  or  a  helix. 
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3.2  GENERAL  ASSUMPTIONS;  INHCMOGENEOUS  WAVE  EQUATION 

To  study  the  propagation  of  waves  in  this  system,  Maxwell's  equations 


VxS--|| 

(3.1) 

VxS-r.ll 

(3.2) 

V*  D  =  p 

(3.3) 

0 

II 

t> 

(3.i^) 

will  be  used.  The  convection  current  density,  1  =  pv  ,  is  assumed  to 
be  due  to  only  the  electrons,  of  the  beam  and/or  of  the  plasma.  It  is 
assumed  that  all  quantities  have  an  average  value  plus  a  small  harmonic 
time  dependent  perturbation: 

F(r  ,  6  ,  z  ,  t)  =  T^{r  ,  9  ,  z)  Fj^(r  ,  0  ,  z)e^“*  ,  (3-5) 

and  the  perturbations  are  wave-like  in  nature,  so  that  Fj^(r  ,  ©  ,  z) 
has  the  form 

P^(r  ,  0  ,  z)  =  Fj^(r  ,  0  )e’'^^^  ,  (3.6) 


where  ^  is  the  propagation  constant  along  the  axis  of  symmetry  of  the 
system.  We  can  write  Maxwell's  equations  for  the  ac  quantities  in  terms 
of  the  spatial  coordinates  as 

Vx  (3.7) 


Vx  =  I 

(3.8) 
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(3.9) 
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where  >  6  >  2)  ,  >  0  >  z)  >  etc.  If  we  now  taJce 

the  curl  of  Eq.  (3.1)  and  substitute  in  Eq.  (3.2)i  we  find  that 

Vx  V  X  , 

where  k  =  ”  •  Using  a  vector  identity  and  Eq.  (3*9) 

Vx  Vx  -  -  v^^*  V(  V-  V 


‘0 

yields  the  Inhomogeneous  wave  equation  for  E^  ,  which  is 

+  jojM  t  .  (3.11) 

*0 


This  equation  can  also  be  written  in  the  fora 


'  •'^^^1 


(3.12) 


r7  2 

where  the  8yiid)ol  v  1  denotes  the  transverse  Laplaclan. 
coordinates,  V 1  »ay  te  written  as 


h  f  h\ 

m  —  —  (X*  —  )  4.  _ 

r  dr  \  dr/  r^  dC^ 


For  cylindrical 


(3.13) 


It  follows  slMlenrly  that 

-  (p2  -  »  0  .  (3.1^^) 

2  2 

If  we  make  the  slow-wave  approxlaatlon  (k  «  0  }  ,  the  dlfferentieil 

equations  satisfied  by  £  and  B  eu:e,  respectively, 
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and 


=  0  . 

z  z 


(3.16) 


The  continuity  equation  is 

V-i‘  +  |^  =  0  .  (3.17) 

In  cylindrical  coordinates,  Eq.  (3*17)  may  be  written  as 


f  (rij.)  -  JPi^  =  -  .'op  .  (3.18) 

This  equation  is  used  to  express  the  volume  charge  density  p  in  terms 
of  the  components  of  the  c\irrent  density,  i.e., 
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(3.19) 


Then  Eq.  (3*19)  nay  be  substituted  into  the  wave  equation  [Eq.  (3*15)]> 
so  that  the  inhomogeneous  term  in  the  equation  becomes  a  function  only  of 
the  current  density: 
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(3.20) 


If  we  restrict  ourselves  to  circularly  symnetrlcal  solutions,  i.e., 
(d/5  6  3  O)  ,  Eq.  (3.20)  becomes 
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This  cylindrical  wave  equation  is  to  be  solved,  in  the  two  regions  of 
Interest  of  the  system  as  shown  in  Fig.  1: 

a.  Region  I  -  (O  <  r  <  a)  is  the  electron  beam-plasma  region, 
defined  to  the  limits  of  the  edge  of  the  beam.  In  this  region  we  find 
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(3.22) 


b.  Region  II  -  (a  <  r  <  b)  is  outside  the  bean  region,  where  only- 
plasma  is  present,  and 

i  =  i 

Zii  zpjj 


(3.23) 


3.3  ELECTRON  BEAM-PLASMA  INTERACTION  REGION 

Consider  the  beam-plasma  region  (Region  I),  (0<r<a)  .  Under 

the  small-signal  assumption,  the  axial  and  radial  components  of  the 
linearized  beam  current  densities  are,  respectively. 
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The  equations  of  motion  for  the  electrons  in  the  beam  lead  to  the 
following  expressions  for  the  beam  velocity  components: 


in 
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where  q  =  e/m  and  the  index  (l)  denotes  the  region  which  is  being 
considered.  From  the  continuity  equation,  applied  to  the  electron  beam, 
we  find  that 
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which,  by  using  the  previovis  equations  for  the  current  densities  and  ac 
velocities,  assumes  the  following  form: 
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v/here 


(3.30) 


is  the  plasma  frequency  of  the  beam.  We  know,  however,  that  with  the 
slow-wave  approximation  the  radial  component  of  the  electric  field,  as 
derived  from  Maxwell's  equation,  is 

E  =  i  ^ 

^  P  Sr 


It  follows  that 


J'o  > 


^  (o)  -  UqP)  r  Sr  '  S 


iSe  \  2 

it)  ■  (0.  • 


By  substituting  the  above  equation  into  Eq.  (3*24)  and  using  (3-26)  we 


obtain  the  result 
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The  radial  component  of  the  beam  current  density  is 


.  0  pb 

^rb  “  ^ObVo  '  ^  / 


2  SE 


P  (o)  -  UqP)  Sr 


(3.32) 


For  the  plasma  electrons  in  the  same  region  (O  <  r  <  a),  the  current 
density  components  are 


(3.33) 
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where 


(3.35) 
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is  the  plasma  frequency  of  the  plasma.  In  this  region  I,  the  total  longi¬ 
tudinal  and  transverse  current  densities  are,  respectively, 


and 


(3.36) 

(3.37) 


Therefore,  by  using  the  expressions  above,  derived  for  the  current  densities, 
we  have  ,  , 


By  substituting  these  last  two  equations  into  the  wave  equation 
we  obtain 


(O)  - 


0  . 


[Eq.  (3.2)], 

(3.IK)) 


We  conclude,  then,  that  there  are  two  possibilities  for  the  fields.  In 
the  beam-plasma  region,  namely: 


(3.1^1) 


(tt 
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(O)  -  UqP)‘ 
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(3.42) 


Consider  now  the  Region  II  (a 
In  this  region,  there  is 
the  current  density  we  write 


<  r  <  b)  . 

only  plasma  and  for  the  two  ccsoponents  of 


1  =1 

"ll  ^PlI 
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(3.43) 


In  this  case  it  may  be  seen  that  the  wave  equation  [Eq.  (3.2l)],  assumes 
the  following  form: 


/  0) 


^  ^  I  2 

- Ir  - —  j  -  p^E 

r  3r  '  Sr  / 


=  0 


(3.H) 


which,  for  cn  o)  ,  is  reduced  to 
’  Pp 


=  0 


II 


(3.J^5) 


3.^  BOUNDARY  CONDITIONS  AT  THE  EDGE  OF  THE  BEAM 

We  calcvilate  the  charge  perturbation  at  the  surface  of  the  electron 
beam,  due  to  the  radial  motions  of  electrons  of  the  beam  as  well  as  of 
the  plasma.  The  ripple  at  the  boundary  of  the  electron  beam- plasma  is 
approximated,  by  the  method  Introduced  by  Hahn,^  by  a  sheet  of  ac  surface 
charge  at  the  same  boundary  the  electron  beam  has  in  the  absence  of  ac 
excitation. 

The  net  surface  charge  density,  at  r  =  a,  results  from  contributions 
due  to  radial  motions  of  electrons  coming  from  region  I  (beam  and  plasma), 
and  also  of  electrons  coming  from  region  II  (plasma  only). 

The  surface  charge  densities  are  calculated  fran  the  small  radial 
displacements  of  electrons  at  r  =  a  .  The  continuity  equation  is 
applied  in  a  small  element  of  volume  on  the  dc  boundary  of  the  beam.  The 
equivalent  surface  charge  density  which  is  due  to  radial  displacements 
of  electrons  of  the  beam  turns  out  to  be 


‘'b  =  pQb 


rb 


J(a>  -  UqP) 
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r  =  a 


(3.i^6) 


r  «  a 


Upon  substituting  in  this  equation  the  expression  derived  before  for  the 
radial  component  of  the  beam  cvirrent  density  (Eq.  3-32),  we  obtain 
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In  a  similar  way,  the  surface  charge  density  at  r  =  a  resulting  from 
plasma  electrons  coming  from  region  I  is  calculated  in  terms  of  the 
electric  field,  yielding  the  eqioation 


=  -  J 


a 


_0  _E2_ 


(3A8) 


Finally,  the  s\irface  charge  density  at  r  =  a  due  to  plasma  electrons 
coming  from  region  II  is 


8e 

"II 


r 


s  a 


(3.49) 


We  know  that  at  r  =  a  the  normal  component  of  the  electric  field, (l.e., 
the  radial  component  of  is  discontinuous  by  the  amount  of  the  total 
equivalent  surface  charge  density  associated  with  the  i>erturbatlon.  We 
have 


(3.50) 


The  other  boundary  condition,  at  r  «  a  ,  states  the  continuity  of  the 
longitudinal  electric  field,  namely 

E*E  ,atr=a.  (3*51) 

"l  "ll 


With  the  help  of  the  expressions  for  the  equivalent  surface  charge  densities, 
as  presented  above,  these  two  bovindary  conditions  lead  to  the  following 
equation: 


hr 
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It  should  be  pointed  out  here  that  the  analysis  made  so  far  applies  each 
set  of  field  solutions  in  the  beam  region,  namely 

a.  the  solutions  deriving  from  Eq.  (3 *^2),  rewritten  below: 


1 


CO 


0) 


.  0 


(co  -  U^^p)^  CO 


and 

b.  solutions  which  result  from 


V 


=  0 


3.5  NONSOLENOIDAL  FLOW  EQUATIONS 

Let  us  consider  now  the  case  (a)  v;here  the  longitudinal  propagation 
constants  p  for  the  fields  in  region  I  are  the  solutions  of  Eq,  (3.^2), 
that  is 

=  B,  i  h  (3.53) 

where  P  =  co/u_  is  the  beam  v/ave  number  and  h  is  defined  as  given  by 
e  o 


(3.54) 


The  pair  of  waves  corresponding  to  these  propagation  constants  is  referred 
to  as  the  nonsolenoldal  mode.  The  electric  field  inside  the  beam  region 
(region  l)  has  nonvanishing  divergence.  Consequently  the  propagation  of 
this  'pair  of  space-charge  weaves  is  associated  with  a  charge  bunching 
(p  ^  o)  within  the  beam.  The  propagation  constants  are  the  same  as  in 
the  case  of  an  infinite  cross-section  beam-plasma,  and  are  therefore 
independent  of  the  finite  geometry  of  the  system. 

Inside  the  electron  beam,  we  have 

Pb^  /  0  3  pp^  0  •  ^3*55) 
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Assume  thc.t  the  modulation  frequency  of  the  beam  <«  is  different 

from  the  plasma  frequency  of  the  plasma  cn  .  (oj  o)  )  .  It  follows 

PP  PP 

then,  from  Eq.  (3.52),  tlmt  the  follov/ing  boundary  condition  exists  at 
r  =  a  : 

=  0  ,  (o)  oSpp)  .  (3.56) 

r  =  a 

We  have,  therefore,  in  the  outside  beam  region  (region  II),  the  longitudinal 
fields  satisfying  the  differential  equation 

eind  the  bovmdary  conditions: 


E  (b)  =  0 


^r 


(3.57) 


The  solution  of  the  differential  equation,  for  a  <  r  <  b  ,  is  written 


E^Il(r)  =  +  B2KQ{^r)  ,  (3.58) 

where  and  are  constants  to  be  determined  with  the  aid  of  the 

prescribed  boundary  conditions.  We  have 

=  p(A2l^(pa)  -  B2K^(pa)j  =  0  (3.59a) 

r  =  a 
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and 


(3.59t>) 


A2lQ(pb)  +  B^K^Cpb)  =  0 


for  each  and  have  to  solve  then,  for  p  /  0  ,  this 

system  of  linear  homogeneous  equations  for  the  unknowns  and  , 

respectively.  It  is  well  known  that  this  system  has  a  nontrivial 
solution  (A^  0  ,  /  O)  if  the  determinant  of  the  coefficients 

is  zero.  This  means  that  we  should  have 


A  =  I^(pa)KQ(pb)  +  lQ(pb)K^(pa)  =  0  .  (3.60) 

However  for  real  values  of  p(p  O),  this  detenninant  does  not  vanish. 

It  follows  then  that  =  0  ,  =  0  .  This  result  implies  that  there 

is  no  field  in  the  plasma  region  outside  the  electron  beam.  The  solution 
in  region  II  is 


(1,2) 


=  0 


(3.61a) 


corresponding  to  the  nonsolenoidal  fields  inside  the  beam.  From  the 
boundary  condition  at  r  =  a  ,  as  expressed  by  Eq.  (3.51),  it  follows 
that 


E  (a)  =  0 
''l 


(3.6lh) 


It  had  been  shown  in  the  preceding  lines  that  when  cu  >  0)^^  ,  these 

nonsolenoidal  modes  do  not  have  fields  external  to  the  beam  region. 

Therefore  they  cannot  couple  to  the  fields  of  an  external  structure. 

However  if  this  structure  penetrates  into  the  beam  by  means  of  an 

obstruction  as  in  the  case  of  grids,  v;hich  will  concern  us,  these  solutions 

become  possible  solutions  for  the  problem.  For  the  case  cu  <  oj  , 

PP  ' 

3  will  have  complex  values;  but  a  similar  consideration  of  (3.60)  for 
this  case  indicates  no  solutions  and  therefore  requires  A  =  B  =  0  and 
Eqs.  (3.6la)  and  (3.6lb)  also  apply. 


-  29  - 


3.6  NONSOLENOIDAL  AND  NONGROWING  WAVE  SOLUTIONS 


Under  the  assumption  that  the  modvilation  frequency  of  the  grid- 
modulated  beam  is  always  above  the  plasma  frequency  of  the  plasma,  we 
have  seen  in  the  previous  section  that: 

a.  There  are  two  possible  propagating  waves  that  can  be  excited, 
one  with  a  phase  velocity  faster  than  ,  and  the  other  with  a  phase 
velocity  slower  than  . 

b.  There  is  no  field  in  the  plasma  region  outside  the  beam.  The 
analysis  which  follows  in  this  section  will  consider  only  this  pair  of 
nonsolenoidal  solutions.  The  treatment  will  be  nonrelativistic,  and  this 
means  that  the  forces  on  the  beam  electrons  caused  by  the  beam  magnetic 
fields  may  be  neglected  and  that  the  electric  field  of  the  electron  beam 
may  be  expressed  as  the  gradient  of  a  scalar  potential  satisfying  Poisson’s 
equation.  This  type  of  description  for  the  space-charge  field  is  some¬ 
times  called  the  quasi-static  approximation. 

Inside  the  beam-plasma  region,  which  has  been  denoted  as  region  I, 
we  have 


■(r,z) 
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(r)e 


-Jhz 
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(3.62) 


where  v^(r,z)  is  the  electric  potential  defined  in  this  region,  and  the 
superscripts  (l)  and  (2)  indicate,  respectively,  the  correspondence  with 
each  nonsolenoidal  pr(^>agatlng  wave  described  in  £q.  (3*4^2). 

The  longitudinal  and  radial  components  of  the  electron  beam  velocity 
are,  respectively. 


and 
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Similar  expressions  hold  for  the  plasma  electrons  velocity  components 
v^^  and  v^  ,  respectively. 
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The  equation  of  motion  for  the  electrons  in  the  besm,  along  the 
z-direction,  ij 


^zb 


-  V 
m  oz  I 


from  which  the  longitudinal  velocity  components  are  obtained  in  terms  of 


the  potential: 


zb^ 


u. 


3e\  V. 
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h  / 
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2V, 


The  radial  velocity  components  are  derived  from  the  equation  of  transverse 
motion,  resulting  in  the  following  expressions: 

^0 

v(2) 

u^ 


(3.64) 


For  the  plasma  electrons,  in  this  same  region  (0<r<a),  we  find  that 
the  longitudinal  velocity  components  are 
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0 


(3.65) 
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and  the  radial  velocity  components  are 


(1) 

rPi 


0 


(3.66) 


Sy  using  the  above  expressions,  and  suppressing,  for  convenience,  the 
phase  factor  ^  we  may  write,  for  the  electron  beam  velocity 

components , 


and 


i 

-  C08(hz) 
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dr 


1 

-  sin(hz) 
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The  corresponding  velocity  components,  for  the  plasma  electrons,  in  the 
same  region  I,  are  given  by 
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-  J  8in(hz) 


-  32  - 


and 


V  J 

— t  =  —  .  cos(hz) 
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It  had  been  shown  in  the  preceding  section  that,  in  region  II 
(a  <  r  <  b),  we  oust  have  Laplace's  equations  satisfied: 


II 


(1)  -y^T 


II 
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(3.67) 


by  the  trivial  solutions: 

=  0  ,  =  0  .  (3.68) 


From  the  continuity  of  the  potentials,  at  the  boxondary  of  the  two 
regions  I  and  II,  the  following  conditions  result,  for  the  potentials,  in 
the  region  I  at  r  »  a  : 

Vj^^^(a)  =  =  0  .  (3.69) 


The  surface  charge  densities,  at  r  >  a  ,  are  now  calculated  for 
the  nonsolenoldal  solutions.  The  contribution  to  the  total  surface 
charge  density,  at  r  =  a  ,  due  to  the  radial  motion  of  the  electrons 
of  the  beam  is  given  by 
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(3.70) 
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where,  for  simplicity,  we  again  omitted  the  phase-factor  e 
We  have 
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a  (2)  V  (2) 
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dr  \  2V, 


(3.71) 


expressions  which  follow  from  equations  (3- 6^). 

The  contribution  to  the  total  ac  surface  charge  density,  at  r  =  a  , 
which  is  due  to  the  plasma  electrons  caning  from  region  I,  is  written 
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-Jhz 
e  "  + 
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where 


‘  \2Vq 


r  »  a 
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3.7  SOLENOIDAL  FLOW;  APPROXIMATE  SOLUTI(»S 

For  the  sedond  possibility  given  by  (3.*+0)  we  assume,  alternatively, 
that  Eq.  (3.^1)  Is  sdtiAfied  in  the  beam-plasma  region.  The  field  solu¬ 
tions  are  therefore  to  be  derived  from  the  Laplacian  of  the  longitudinal 


electric  field 


E  =  0  . 


(3.73) 


If,  however,  we  fefer  to  Eq.  (3.21)  and  write  the  right-hand  side  of  this 
equation  in  terms  of  the  charge  density  with  the  help  of  the  continuity 
equation,  we  have 
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It  follows  then,  from  Eq.  (3 •73),  that  we  must  have 

=  0,  (p,  =  0  ,  p  =0)  ,  (3-75) 

I 

and,  consequently 

div  .  =  0  .  (3-76) 

We  refer  to  these  modes  as  solenoidal  modes  in  contrast  to  the  other 
solutions  considered  before  -  nonsolenoidal  solutions  -  where  the  fields 
have  nonvanishing  divergence. 

In  this  solenoidal  case,  where  the  divergence  of  the  electric  field 
is  zero  everywhere  within  region  I,  there  is  no  ac  voliane  charge  density 
inside  the  beam.  These  solutions  correspond  to  space-charge  wave  pro¬ 
pagation  involving  a  rippling  at  the  surface  of  the  electron  beam. 

The  fields  produced  by  these  modes  are  caused  entirely  by  the  bulging 
and  contracting  of  the  surface  8ex>aratlng  regions  I  and  II. 

In  solving  the  differential  equation 

-  (0  <  r  <  a) 

we  assume  solutions  of  the  form 

E^^(r)  =  Aj^lQ(3r)  ,  (3.78) 

where  a  second  solution,  which  is  not  finite  at  the  origin,  has  been 
omitted. 

In  region  II  (a  <  r  <  b)  ,  we  also  have 

=0  ,  (3.79) 

V  Zjj 

and  the  solution  of  this  equation  satisfying  the  prescribed  boundary 
condition  at  r  =  b  ,  is 

®Zjj(r)  =  A2(lQ(pr)KQOb)  -  lQ(pb)KQ(pr)]  ,  (3.8o) 

for  a  <  r  <  b  . 
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By  using  Eq.  (3*52)  with  the  above  expressions  for  the  fields  in  the 
two  regions,  the  following  determinantal  equation  is  obtained: 
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0)  ,  0) 
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I^(Ba)KQ(pb)  + 
lQ(pa)KQ(pb)  -  lQ(pb)KQ(pa) 


from  which  we  get 


Il(pa)  lQ(pb)KQ(pa)  -  Ip(pa)KQ(pb) 
iQ(pb)  *  I^(pa)KQ(pa)  +  ypa)K^(pa) 


>  .(3.81) 


With  the  aid  of  the  WIronsklan  relation, 


l3^(pa)KQ(pa)  +  lQ(pa)Kj^(pe)  -  —  , 


we  write 
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(3.82) 


By  making  the  identifications 
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we  may  write 


p  =  Pg  +  hr(p)  .  (3.85) 

This  detenninantal  equation  is  a  rather  involved  equation  to  be  solved 
for  p  .  Instead  of  trying  to  solve  it  in  the  most  general  case,  we 
proceed  by  making  simplifying  assumptions,  which  are  shown  below  to  be 
valid,  under  a  practical  point  of  view. 

To  begin  with,  we  assume  that  o)  >  0)^^  ,  which  means  that  we  shall 
be  concerned  with  real  values  for  p  . 

We  denote  the  propagation  constants  for  this  pair  of  solenoidal  waves 
by 

P3  =  pg  +  hF3(P3)  =  Pg  +  A3  (3.86) 


and 

PU  *  Pe  -  =  Pe  ■  ^  ’ 


where,  by  definition 


A3=  hF3  =  hF(p3) 

and 

A4  =  hFj^  =  hFO]^)  • 

We  assume  further  that 

^<1  .  (3.87; 

Pe 

The  degree  to  which  such  approximation  can  be  applied  can  be  estimated 
by  using  some  typical  figures.  In  the  experiment  to  be  described  in  a 
subsequent  chapter,  we  have 
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We  should  have  then,  apart  from  the  factor  P,  the  following  Inequality 


<  1 


0.  2 

_2E_  , 
2  ^ 
cu 


which  leads  to  the  conclusion  that  this  condition  is  Indeed  fulfilled 

for  values  of  o)  very  near  co  (o)  ~  99^  u))  . 
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With  this  approximation,  we  have  for  a  first-order  solution 


(3.88) 


where  A  is  defined  as  being  given  by 


A  .  hPO^) 


(3.89) 


Figure  2  shows  F(^g)  $  which  is  always  less  than  unity,  as  a  function 
of  B'Hd  b/a  .  nils  factor  F(p^)  is  the  usual  "plasma  frequency 

reduction  factor"  which  takes  into  account  effects  of  space-charge  on 
bunching  in  klystrons. The  curves  in  Fig.  2  were  taken  from  Feenberg's 
paper  (ref.  6). 

Under  this  first  order  approximation,  we  are  allowed  to  write 


+  hF(p^)  -  p^  +A 

p,^  -  P^  -  hF(p^)  -  P^  -A  .  (3.90) 

As  a  next  approximation,  we  make 

p^  «  P^  +  hP(p^  ♦a)  -  P^  +  A' 

p,^  -  p^  -  hF(p^  -A)  -  P^  -A"  .  (3.91) 

It  follows  that 
A'  «  hF(p^  +  a) 


1  ^ 

.  h  ^F(P^)  +  A  -  .  F(p^)|>  =  hF(p^)  +  h^F(p^)  .  - F(p^)  , 
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Fig,  2 — Plot  of  the  reduction  factor  F  for  different  values 
of  p^a  (taken  from  Feenberg's  work). 
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or 


A'  -  A.h%'(S  )  .  —  F(f,)  , 

®  dp  * 


and  similarly 


We  have,  then. 


+  A'  »  t  hr(p^)  th®  r(e^)  .  —  F(p  )  , 


and 


»  P^  -  hF(p  )  +  h^(p  )  .  - F(6  )  , 


from  which  we  identify 


and 
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By  putting 


F^(P^)  =  G(pg)  , 

we  have 

6  A-  -  A,,  Sg 

^ - - .  (3.96) 

ha  ha  ^Cp^a) 

Figure  3  shows  SG/5(p^a)  as  a  function  of  (Pg®')  b/a  . 

These  curves  were  adapted  from  the  curves  presented  by  Feenberg. 

3.8  MATCHING  INPUT  BOUNDARY  CONDITIONS  WITH  THE  SPACE-CHARGE  WAVES 

In  matching  the  input  conditions  with  the  space-charge  waves  discussed 
in  (3.^  and (3. it  is  convenient  to  express  the  fields  for  the  solenoidal 
modes  as  the  gradients  of  potentials  Vj  (r,t)  and  (r,  z)  that  ful¬ 

fill  Laplace's  equations,  respectively, . inside  and  outside  of  the  electron 

beam.  This  is  valid  for  nonrelativistic  velocities. 

2 

To  determine  the  solutions  of  y  =  0  in  region  I,  as  the 
solutions  appropriate  to  the  solenoidal  case,  where  the  propagation 
constants  have  the  approximate  expressions  derived  before,  we 

assume 

Vj<3)(r)  =  Ajypjr) 

Outside  the  beam  reglorii  we  have 

^2  V  =  0  >  (3.98) 

V  II 

the  solutions  of  which  must  satisfy  the  following  bovindary  conditions: 

=  0  (3.99) 


-41- 


and 


.  (3.100) 

The  solutions.  In  the  region  II,  then,  are 

=  B^[lQ(p^r)KQ(3^b)  -  lQ(p^b)KQ(p^r)] 

=  Bi^[lQ(p^r)  .  yp^b)  -  lQ(a^b)y P^r) ]  ,  (3.101) 

vhere  the  constants  and  Bj^  must  be  det  rained  in  order  that  the 

prescribed  boundary  conditions  at  r  =  a  ,  are  fulfilled.  With  the  help 
of  these  conditions,  and  Eqs.  (3.100),  ve  obtain 

Aj  .  (Io(»3.)K„(S3b)  -  I^OjtjKodJy.)] 

and 

\ 

\  - - - -  tlQ(pj,a)  .  KQ(Pj^b)  -  lQ(Pi^b)KQ(pj^a)]  .  (3.102) 

Additional  relations  among  the  constants  arise  from  one  of  the 
conditions  at  the  input  plane  (z  =  O).  We  should  have  ^];j(^)  ~  const, 
s  0  ,  for  &  5  z*  <  B  .  This  imposes  the  condition  that 

Vjj(a)  =  0  .  (3.103) 


That  is. 


'II 


^^^(a)  +  Vjj^^^(a)  =  0 


(3.103b) 
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since  the  nonsolenoldal  fields  within  the  beam  do  not  have  fields  in 
region  II.  This  equation  together  with  Eqs.  (3*97)  yields  the  following 
relation  between  the  coefficients  and  Aj^  : 


=  -  Aj^I^Cpi^a)  . 

Note  that  condition  (3*103)  results  in  *  0  for  a  <  r  <  b  only 

within  the  approximation  that  3^  «»  or  3^/h  »  1  .  The  solenoldal 

potentials  may  then  be  written  in  terms  of  one  of  the  coefficients  A^ 
or  Aj^  : 


Io03a) 


=  A^lQ(3^a) 


(3*104) 


The  nonsolenoldal  field  solutions  which  have  the  same  propagation  constants 
as  in  an  infinite  system  are  allowed  to  have  an  arbitrary  transverse 
dependence. 

We  proceed  by  applying  the  remaining  input  boundary  conditions  to 
the  two  pair  of  waves  which  are  excited  in  the  grldded  gap.  The  radial 
variation  for  the  solenoldal  modes,  as  shown  above,  is  completely  specified 
from  the  solution  of  the  Bessel  differential  equation  [Eq.  (3*73)1  inside 
the  beam. 

The  boundary  conditions  to  be  satisfied  at  the  input  plane  will 
impose  then  a  particular  transverse  dependence  on  the  nonsolenoldal  modes. 
This  dependence  is  not  the  most  general  radial  vsu:iation  which  could  be 
prescribed  for  these  nonsolenoldal  field  components.  Feenberg  shows  that 
this  general  solution  includes  terms  which  are  negligible  under  spproxi> 
matlons  previously  made  in  deriving  the  propagation  consteuits  [terms  of 
order  and  higher  order]. 

However,  all  that  is  really  required  is  to  show  that  the  radial 
variation  used  here  for  these  waves,  together  with  the  solenoldal  waves, 
fulfill  the  boundary  conditions  of  the  problem. 


In  the  remainder  of  this  section  this  relation  between  the  two  pair 
of  modes  is  determined. 

At  the  input  plane  (z  =  O)  ,  by  assuming  ideal  grids,  we  must  have 
the  following  condition  inside  the  beam  region: 

Vj.(r)  =  0  ,  at  z  =  0  . 


The  input  longitudinal  velocity  bf  the  electron  beam  is  a  specified 
quantity.  Let  v^^  *  ^0  z  =  0  .  Hierefore  we  may  write 


il 


(3.107) 


where  is  the  dc  beam  velocity,  M  *  2)! 2)  beam 

coupling  coefficient,  and  a  is  the  voltage  modulation  index  (a  » 

By  using  Eqs.  (3.63)  and  the  similar  ones  for  the  solenoldal 
(3)  (4) 

components  ^^b  \b  beam  longitudlnail  velocity, 

we  obtain  the  following  equation: 


This  equation  together  with  Eq.  (3.106)  yields 


QhM  IpOi^a)  /  1  Io(P4r)  1 

2Vo  “  "  2p/  2V^,  *  ^  F3  '  ' 

(3.108) 

But,  from  the  continuity  of  E  at  r  s  a  and  the  fact  that  E  *  0 

%  z 

for  r  >  a  for  the  nonsolenoldal  waves. 


Vj^^^a)  -  Vj^2\a)  =  0  , 


and  It  follows  that 


» 


and 


(3.109) 


We  have,  then,  the  solenoldal  components  of  the  electric  potential  In 
the  beam  region  given  by  the  following  expressions: 


io(»3») 


and 


Vo^V^ 


W) 

Io(P4a) 


(3.110) 
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The  other  input  boundary  condition  to  be  satisfied  is 


v,=0  ,  at  z  =  0  . 

rb  ’ 


^3)  (l^,) 

Using  Eqs.  (3-64)  and  similar  ones  for  '  ,  the  following 

equation  results: 


d  j  1  d 


1  d  /V 


^3  \  2Vo  /  F4  dr  \  2Vq 


=  0  , 


which  can  be  written 


d 


QM/  h 


I  (p-r)  Pr 
-± — 2 —  +  - 


V.,..  ,(,.L  vv. 


where  we  have  used  Eqs.  (3.109)« 

The  condition  on  the  longitudinal  current  density  of  the  electron 
beam  at  the  entrance  plane  z  =  0  is  stated  as  follows: 


(3.111) 


We  know,  however,  that 


(1) 

zb 

'  ''V 

(2) 

zb 

•  “oj^fob  ■ 

,  v(2) 

zb 

(3) 

s  ^ 

„(3) 

zb 

^Ob  * 

^zb 

w 

v('^) 

zb 

“  ^Ob  • 

^zb  ' 

(3.112) 
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therefore, 


.(1)  ..(2) 


(3.113) 


Moreover,  the  volume  charge  densities  associated  with  the  nonsolenoidal 
fields  can  be  calculated,  in  terms  of  the  corresponding  potentials,  by 
applying  the  continuity  equation  at  z  =0.  We  have 

d 

—  i  ,  =  -  Jop.  ,  at  z  =  0  . 


Then  we  find 


O'  *  ‘zi’  ■ 

■  O'  ■  0.)  ‘il’  •  “Ob^’  • 

5y  using  Eqs.  (3.112)  for  i^^^  and  i^^^  ,  and  substituting  in  the 
above  equations,  we  have 

p(l)  / 

% - 1  .  ^  ^ 

\  h/ 


making  use  of  Eqs.  (3*63}  ve  can  express  '  and  '  in 

(1)  (2)  ®  ^ 

terms  of  the  potentials  Vj  '  and  p  respectively,  yielding 

(1)  /  ^  \2  M) 


-  l'^- 

\  h 


-  k8  - 


and 


(3.11^) 


3.9  APPROXIMATE  SOLUTIONS  FOR  THE  SPACE- CHARGE  WAVES;  VOLUME  AND  SURFACE 
CURRENTS 

In  this  section  we  restrict  ourselves  to  the  approximations  for  the 
propagation  constants  of  the  space-charge  wave  modes,  as  indicated 
in  Section  3* 7* 

We  assume  then  that  the  following  inequalities  remain  valid: 

a.  CD  >  (u 
PP 


Under  these  approximations,  we  have 


(3.115) 


±  hP(p^)  =  P^  ±  A  ,  13.116) 

and  we  can  expand  the  functions  and  lQ(Pj^r)/lQ(Pj|j^a)  , 

respectively,  in  power  series  in  A  ,  to  obtain 
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and 


Substituting  these  two  expansions  In  Eqs.  (3*106)  and  (3*108),  respectively, 
results  in  the  follovd.ng  expressions: 


OM  /  A  \  ^ 


2V^ 


2\Q 


L 


+  0  - 

\Kj 


and 


M  y(2) 

I  *I 
2V« 


CM  ^h 
2 


& 


X.  W, 


ys.o) 


Putting  these  expressions  Into  £q.  (3«l4)  yields 


"Ob 


!\b/  Io(S,a)  \fj 


and,  as  It  has  been  shown  before. 


(3.117) 


(3.118) 


(3.119) 


» pi**' 


"Ob 


CM 

2 


The  nonsolenoldal  cosiponents  of  the  ac  bean  velocity  are,  then,  calculated, 
a.  The  longitudinal  velocity  components  are 

,(1)  .  J2) 


and 

v^l)  -  v<2) 
zb  ^zb 


V'  '  +  v\ 
zb  zb 


u. 


OM  /S  \/A 


CM 

2 

5 

'  .< 


lo(Pea) 


P- 


^  /W^ 

'  ^Pe  VVPe®)^^'^  °ip- 


>  . 


(3.120) 
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b.  Tlie  padlal  components  are 


,(1)  ^  ,(2) 
rb  ^  ^rb 


u. 


QM 
2J 


and 

’rt’  -  'rt' 


-■  ~  +  0  -5-1 

Io(Pe“) 


—  —  < 


U, 


2J\  b  /ve 


3  .  —  '  ' +  0  — 


(3.121) 


In  this  same  region,  the  nonsolenoldal  components  of  the  plasma 
velocity  are  calculated,  yielding  the  following  expressions: 

^(1)  ^  ^(2) 


A 

^0 

2Pe  V 

v(l) 

^Pj 

.  v(2) 

zPi 

OM  /  A\  ^ 

"0 

2  VpJ 

v(l) 

rPi 

^^(2) 

rPi 

OM  / 

"0 

Uey 

,(1) 

.v(2) 

rPi 

OM  /  h  \ 

■-(rj 

"0 

d  Ai(V) 

p  —  A-S-Uo 


cy 


r (3.122) 


Appendix  B  contains  the  derivation  of  the  total  volume  current.  The 
results  are  presented  below.  The  electron  beeun  volume  current,  at  z  *  L 
(output  gap),  is  given  by 


vb  _ 


OM^ 

21,0  a) 

-11-  ( 

J  sin(hL)  + 

e  ^  .  co8(hL)  -  cos(hL) 

2 

\h 

- 

(3.123) 
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Figure  (1^)  is  a  plot  of  21,  (p  a)/(3  a)l-(p  a)  versus  (3  a) 

^  X  ©  0  U  ©  6 

taken  from  Feenberg’s  work. 

The  plasma  volume  current  at  a  distance  z  from  the  input  gap  is 
given  by  the  following  expression: 


0“^  »  /  21,  (pa)  \  /A 

- pE]  J_  [i - ?  .  .co.(ha)  -  J(- 


-j|z  2I3_(^ga) 


sln(Az)  +. . . 


(3.124) 


It  may  be  seen  from  the  above  expression  that  the  amplitude  of  the  volume 
plasma  current  is  considerably  larger  than  the  beam  current,  due  to  the 
presence  of  the  factor  »  1  .  The  surface  currents  are 

calculated  in  Appendix  C.  The  total  surface  current  at  z  »  L  (output 
gap)  turns  out  to  be  given  by  the  expression 


'.b 


aiiO.a) 


Ip  2  '  (P,b)Ip(P,a) 


sln(^)  sln(hL) 


i  ^  COS6L  >  +.. 


(3.125) 

Finally  the  total  current,  at  the  output  grids  (z  =  L),  will  be  given  by 


I  lyb  I.b  -J|l. 

—  - - ^  jp  e  -  .  - 

Iq  *0  ^0  2  *  (Vlio'S.*) 


sin(hL)  21,  (pa)  f -J& 

-  1  -  -  +  e  cos(/^)  -  cos(hL) 

b  (p^a)lQ(p^a) 


2Il(Pe^) 

(p^a)IoOea) 


(3.126) 
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Fig.  4- -Coefficients  which  appear  in  the  expression  for  the  rf 
current,  as  a  function  of  (P^a)  (Feenberg). 
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This  equation  constitutes  the  formal  solution  to  the  problem  of  finding 
the  driving  current  at  the  output  grids.  For  large  values  of  p^a(p^a  >  «) 
this  equation  is  reduced  to  the  expression  for  the  total  current  corres¬ 
ponding  to  the  infinite  cross-section  plasma-electron  beam,  which  is 


I 


QM^  sln(hL)  j(ajt-3  z) 

-  •  JPa  - « 

2  ®  h 


(3.127) 


As  pointed  out  by  Feenberg,  it  may  be  seen,  from  the  expressions  for  the 
volume  and  surface  currents,  respectively,  and  with  the  aid  of  Figs. 

(2)  and  (4),  that  the  first  term  in  Eq.  (3.126)  may  exceed  the  second  by 
a  factor  of  ten.  The  total  surface  current  may  then  contribute  significantly 
to  the  total  current  in  the  output  gap.  However,  it  will  'be  shown  in  the 
following  sections  that,  when  <0  <  0)^^  and  growing  waves  are  then  present 
in  the  system,  the  rates  of  growth  of  the  waves --having  ^olusie  current 
(nonsolenhldal)— are  larger  than  the. growth  rates  of  the  solenoldal  waves. 

3.10  GROWING  WAVE  CASE 


For  the  case  in  which  the  beam  modulation  frequency  is  below  the 

plasma  frequency  of  the  plasma,  it  was  shown  that  the  propagation  constants 

become  conqplex,  indicating  a  growth  in  amplitude  for  the  space-charge  waves. 

We  have,  for  oi  <  o  , 

’  PP  ’ 


and 


(3.128) 
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If  we  nf^nln  oonfjne  ourr-elver.  to  the  range  of  validity  of  the  In- 
'■‘qunl.  U  lofj 


A* 

—  <  1 


the  analysis  mode  in  the  preceding  section  holds.  The  potential  solutions 
in  region  II  outside  the  beam  are  then  given  only  by  the  solenoldal 
components ; 


^(r,z)  .  ^jj^(r)e"'^’’‘  +  vJJ^(r)  •  e*^^®  ^ 


where,  from  Eqs.  (3. 101)  and  ( 3.109),  we  have 


/,\  OfhTW^ 

v[3)(r)  -.j - 0__ 

/I  1 


/LX  ahW. 

vJJ)(r)  -  J  — - ^ 

/ 1  1 


.  Vp3r)Ko(P,b)  -  lQ(P3b)KQ(p^r) 

-  VP3^)Ko03a) 


W)W)  -  Vp4^)W) 

lo(v)VV)  - 


(3.129) 


The  arguments  of  the  modified  Bessel  functions  in  the  above  ejqjtreaslons 
are  complex  quantities. 

The  particular  case  for  which  (u  ■>  corresponds  to  a  rather 
pectillar  sltuatl(xi  at  idiich  the  above  analysis  does  not  hold.  The  fields 
in  the  plasma  region  outside  the  beam  do  not  satisfy  Eq.  (3.45).  This 
means  that  the  expressions  for  the  electric  potentials  in  this  region  as 
given  by  Eqs.  (3.I29)  are  no  longer  true,  since  they  are  not  derivable 
from  Laplace's  equation. 

Moreover,  at  o>  ■  0)^^  the  condition  for  the  fields  on  the  boundaries 
of  the  electron  beam,  as  expressed  by  Eq.  (3-52),  breaks  down;  the  matching 
problem  becomes  unspecified. 
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These  objections  may  be  somewhat  academic  since  the  whole  theory  is 
based  on  a  mathematical  model  of  surface  currents  which  would  not  apply  if 
infinities  occur  in  the  equation.  More  serious  however  is  the  situation 
for  0)  near  0)^^  when  quantities  become  large  even  though  not  Infinite 
and  the  approximation  h/p^  <  1  seems  invalid.  In  this  case,  which  is 
discussed  in  more  detail  in  the  next  chapter,  it  can  be  shown  that  ther^nal 
velocities  in  the  plasma  modify  the  theory  in  such  a  way  that  th«  approxima^- 
tion  VPg  ^  ^  still  valid. 

3.11  CONCLUSIONS 


It  is  of  value  to  summarize  the  results  of  this  chapter  for  the 
finite  electron  beam.  For  the  general  case  of  an  electron  beam  going 
through  a  plasma  with  no  magnetic  field,  the  space-charge  waves  are  of 
two  kinds.  One  set  is  called  the  solenoidal;  it  has  no  space-charge 
density  within  the  /olume  of  the  beam  but  does  have  a  rippled  beam  surface 
which  constitutes  a  surface  current.  The  other  set  is  called  nonsolenoidal; 
these  solutions  have  both  a  volume  charge  density  and  a  surface  ripple. 

The  second  set  is  characterized  by  having  no  electric  field  exterior  to 
the  electron  beam,  aiid  the  charge  distribution  in  the  interior  may  have  an 
arbitrary  dependence  on  the  transverse  co-ordinates.  All  these  results 
are  entirely  analogous  to  results  obtained  for  an  electron  beam  in  the 
absence  of  a  plasma  (and  no  magnetic  field).  The  major  difference,  if  a 
plasma  is  present,  is  that  all  the  propagation  constants  are  modified  by 
the  equi/alent  dielectric  constant  representing  the  characteristics  of  the 
plasma.  The  most  significant  result  is  that,  at  frequencies  below  0^^  , 
the  si>ace-charge  waves  which  can  exist  are  growing  waves.  This  is  true  of 
both  the  solenoidal  and  nonsolenoidal  waves.  With  grid  modulation,  both 
kinds  of  waves  are  excited,  will  exist  in  the  beam,  and  will  contribute  to 
the  current  in  the  electron  beam  at  the  output  gap.  This  is  contrary  to 
the  case  of  a  gridless  gap  where  only  the  solenoidal  solutions  are  produced. 
Another  significant  result  of  this  calculation  is  that  the  growing  waves 
excited  by  grid  modulation  which  are  nonsolenoidal  have  a  greater  growth 
rate  than  the  solenoidal  waves  and  therefore  in  any  amplification  process 
will  predominate. 
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CHAPTER  IV 


HARMONIC  GENERATION  IN  A  VELOCITY -MODULATED  ELECTRON  BEAM 
IN  A  PLASMA  MEDIUM 

k. l  INTRODUCTION 

In  this  chapter  we  shall  use  the  Eulerian  approach  to  derive  the 
nonlinear  space-charge  wave  equation  for  the  second  harmonic  field  com¬ 
ponent.  The  solution  of  this  nonlinear  differential  equation  is  obtained 
and  conditions  for  second  harmonic  growing  waves  are  investigated.  Gen¬ 
eral  expressions  for  the  second  harmonic  quantities  are  derived,  and  some 

of  the  implications  of  the  resultant  formxjdas  are  discussed. 

12 

Paschke  developed  a  nonlinear  space-charge  wave  theory  for  the 
special  case  of  a  drifting  velocity-modulated  beam.  This  theory  was 
approached  by  successive  perturbations,  beginning  with  the  linearized 
solution. 

The  purpose  of  this  chapter  is,  by  using  Paschke *s  technique,  to 
derive  solutions  for  the  infinite  beam  drifting  in  an  Infinite  cross- 
section  plasma.  Unfortunately,  this  method  of  successive  approximations 
became  extremely  lengthy,  for  higher  harmonics,  so  that  we  shall  restrict 
ourselves  to  second  harmonics. 

The  analysis  which  follows  is  based  essentially  on  the  following 
assumptions : 

a.  An  input  signal  of  low  power  modulates  the  beam.  The  beam  electron 
velocities  as  well  as  the  plasma  electron  velocities  remain  single-valued, 

l. e.,  no  overtaking  of  beam  or  plasma  electrons  occurs. 

b.  Infinite  magnetic  field  is  applied  along  the  direction  of  motion  of 
the  electrons,  which  is  taken  as  the  z-axis.  Therefore,  there  is  no 
transverse  motion  of  electrons. 

c.  The  velocity -modulated  electron  beam  is  of  infinite  extent  and  has 
a  dc  velocity  along  z  . 

d.  The  plasma  ions  are  stationary,  i.e,,  are  of  infinite  mass.  Thermal 
motions  and  collision  effects  are  neglected. 

We  assume  also  that  all  ac  quantities  vary  as 

nj(a)t-P  z) 

^z,t)  .  ^  «l„(V)e  *  ,  (4.1) 

n=l 
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where  P  ^  is  the  beam  electronic  wave  nvunber  and  h  denotes  the 

®  ®  th 

electron  beam-plasma  wave  number  characteristic  of  the  n  harmonic. 

We  shall  consider,  in  the  above  equation,  the  case  where  the  funda¬ 
mental  and  the  second  harmonic  components  (n  =  1,2  ,  respectively)  are 
the  relevant  ones,  and  thus  neglect  the  contributions  from  the  higher 
harmonics.  We  may,  therefore,  write 

J(a)t-P  z)  2j(a)t-P  z) 

^(z,t)  =  ^^(h^z)e  +  , 

where  h^^  *  ‘‘^pb^^Opl  ^1  *  ^  defined  for  the  first 

order  variables  (n  «  l)  ,  and 


(U 


oi  2 


(2a))‘ 


(i^.2) 


are  related  to  the  second  order  components  (n  «  2)  .  The  corresponding 
complex  conjugate  of  ^(z,  t)  is  written 


-j(a)t-P  z) 

^(z,t)  -  ^(h^z)e 


^(h2z)e 


-2j(a)t-P^z) 


(i^.3) 


and  we  shall  denote,  for  convenience,  the  svun  and  the  difference  of  ^ 
with  its  complex  conjugate  ^  ,  respectively,  by  the  superscripts  (+) 
and  (-)  . 

4.2  LINEARIZED  THEORY;  FUNDAI4ENTAL  COMPONENTS 

By  means  of  Eq.  (4.1),  the  first  order  ac  variables  are  first  of 
all  defined.  The  second  order  veuriables  are  assumed  negligible.  Linear¬ 
ized  equations  will  result  and  the  well  Imown  results  of  the  linear  space- 
charge  wave  theory  eu*e  obtained.  For  the  sake  of  conqpleteness,  we  shall 
derive  in  this  section  the  second  order  differential  equation  for  the 
fundamenteiL  component  of  the  electric  field.  In  the  following  section 
the  second  haraonic  components  are  considered: 
a.  Continuity  equation  applied  to  the  electron  beam: 


d 

dz 


(4.4) 
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where  the  linearized  beam  current  density  is  defined  as 


‘bl  '  “o\  * 


Ct.s) 


b.  Force  equation  applied  to  the  electron  beam: 


S  S  X 

—  v,  +  u„  —  V.  =  —  D,  , 

^  1  ^z  ^  ^Ob 


(4.6) 


where 


_.-f  _+ 

■>1  ■  '0®1  • 

c.  The  corresponding  equations  for  the  plasma  electrons  eure,  respectively, 


-  i  - - P 

bz  ^1  St  ^1 


(4.7) 


where  the  linearized  plasma  current  density  is  defined  as 


“  ^Op^'p^  ’ 


(4.8) 


S  o) 

_  v  +  .  -EE> 


(4.9) 


Poisson's  equation  reads 


TN+  +  + 


(4.10) 


The  velocity  of  the  plasma  electrons  is  obtained  from  Eq.  (4.9)  as 


cu  2 
-  J 


(4.11) 
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The  plasma  current  density  is,  therefore,  by  using  Eq.  (4.8), 


2 

^  D-  .  (4.12) 

Pi  m  ^ 

With  Eqs.  (4.7)  and  (4.10)  combined  together,  the  result  is 


-  1^1+  —  =  —  Pv  > 

8t8z  8z  Pi  8t  X 

from  which,  after  replacing  8/dt  by  Ju)  ,  we  obtain  the  following 
expression  for  the  electron  beam  charge  density  in  terms  of  the  electric 
field: 


(4.13) 


The  plasma  charge 
equation  for 

Pp  *  (1  -  «i)  —  •  (4.14) 

*^1  dz 


density  is  found  from  Poisson's  equation  and  the  above 
to  be 


The  continuity  equation  for  the  electron  beam,  together  with  Eq.  (4.13), 
yields 


d 

dz 


'^0*1 


Ob 


> 


which  is  substituted  into  the  equation  of  motion,  to  give  the  beam  ac 
velocity  in  terms  of  the  field: 


€  U  S  €  U  d 

^10  -  10  4* 

“Pob  Pob 


2 

-2^  n- 


cup^ 


Ob 


(4.15) 
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With  the  help  of  Eqs.  (4.13)  and  (4.15),  "the  electron  beam  current  density 
is  seen  to  be 


+  ^1^0 

i  =  -  1  - 

b, 

1  O) 


~"o  D  +  2u^e.  —  -  j  '  D. 

^.2  1  0  1  1  1 


°  1  SZ  -  0) 


(4.16) 


It  is  known,  however,  that  in  this  one-dimensional  problem  the  toted 
current  (convection  plus  displacement  current)  must  be  zero: 


(4.17) 


It  follovs,  then,  from  Eqs.  (4.l6),  (4.17)  and  (4.12)  that 


2 

bz^  Uq  BtSz 


d;  *  -  ?/)D- 


(‘*.18) 


The  solution  of  the  homogeneous  differential  equation  corresponding  to 
the  case  of  an  infinite  beam  in  a  drift  space  has  been  derived  by  several 
authors  (see,  for  example,  A.  H.  W.  Beck  ^).  For  the  electron  beam  in 
a  plasma,  the  solution  of  the  differential  equation  (4.l8)  vhich  corres¬ 
ponds  to  the  pair  of  space-charge  waves  with  propagation  constants 

P  *  P  ±  h  ,  is 
e  1 


2, 

^O'^O  ^1 


_2b 

OJ 


cos  h^z  + 


Ur. 


sin  h^z  >  e 


(4.19) 


j(ajt-pgZ) 


where  q  =  e/m  ,  and  i^  (O)  and  (O)  are,  respectively,  the  beam 
current  density  and  the  Hieam  velocity,  at  the  input  gridded  gap  (z  *  O) 


-  61  - 


The  electron  beam  velocity  and  current  density  may  be  shown  to  be, 


respectively, 


(0) 


t, 

—2^ 


k  (0)^ 


l\ 


\  -  •  cos  h^z  -  J 

^  I  ^0  \  c 


j(o)t-P^z) 


sin  h^z  ;>  e 


,  (^.20) 


iv  (0) 


•  cos  hj^z  -  J 


ttj/cT  1  j(a)t-P  z) 

sin  h,z  S  e 


O) 


(^.21) 


-0  pb 

The  boundary  conditions  to  be  satisfied  at  the  input  gap  are 


S^(0)  .  0 

v^^(O)  “  *  UqO  I  ,  (4.22) 

D^(0)  .  0 

vhere  a  «  depth  of  modxilation,  and  M  denotes  the  gap 

coupling  coefficient. 

Under  these  boundeiry  conditions  the  expression  for  is  reduced 

to 


j(a)t-B  z) 


sin  h^z  •  e 


I  i 


from  which  we  have 


and 


21  /an\ 

D^(z,t)  «  +  1^  - - I  —  j  •  sin  hj^z  •  cos  (cut-B  z) 

“ptfiVj 

‘o 

D,(z,t)  *  D,  -  D*  *  2J  - I  —  1  ■  sin  h.z  •  sln(cot-B  z) 


(4.23) 
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The  electron  beam  velocity  and  current  density  become^  respectively, 

QM  j((ot-p  z) 

V,  =  u-  —  •  cos  h,z  •  e  ®  (4.24) 

bi  0  2  1 


and 


4.3  FIELD  SOLUTIONS  FOR  THE  SECOND 
SOLUTIONS 


j(<nt-6  z) 

sin  h^z  •  e  .  (^.25) 

HARMONIC  COMPONENTS;  GROWING  WAVE 


In  this  section,  we  shall  Include  the  second  order  terms,  in  the  form 
of  Eq.  (^.l),  in  the  expressions  for  the  ac  quantities. 

These  second  order  variables  are  solved  in  terms  of  the  fundamental 
components  which  were  determined  in  the  previous  section. 

We  start  by  defining  the  second  harmonic  components  of  the  electron 
beam  and  plasma  current  densities,  which  are,  respectively, 


.+  +.  +.lt  + 

^bg  “  Vb2  P0b^b2  2  Pb^^b^  ' 


(4.26) 


and 


.+  +  1  +  + 
i  -  p-  v  +  o  P  V 
Po  Op  Po  2  p^ 


(4.27) 


^2  *'2  “  ■*'1  •‘'1 
For  the  electron  beam,  the  continuity  and  force  equations  are,  respectively. 


St 


(4.28) 


and 


—  o  ■''’v  '  —  V, 

ot  ^2  °  Sz  ^2  ^  ^1  Bz  \ 


Ob 


The  corresponding  equations  for  the  plasma  electrons  are 


(4.29) 


a 

Sz 


(4.30) 


-  63  - 


and 


d  -  S  0) 

+  ,  1  +  +  pp  _+ 

bt  ^2  ^  bz  ^1  p  ^ 

Op 

Poisson's  Equation  applied  for  the  second  harmonic  components  Is 

^  n*  *  * 

-  Do  =  Pv  + 


(^.31) 


(^.32) 


Take  the  ptirtlal  time  derivative  of  the  above  equation  and  use  Eq.  (4.30) 
to  obtain 


at  ^^2 


-■  I  -’>S*  Pop\  *  2  W 


dz 


ht 


(4.33) 


/ 

After  eliminating  v  from  Eqs.  (4.31)  and  (4.33)  we  find  that 


€p  •  —  Dp  -  — ^  •  — 

dz  ^  dz 


1  «+  + 

dt{^ 


^Op  +  ^  + 

-  V  —  V 

2  ^1  bz  H 


(4.34) 


The  continuity  equation  and  equation  of  motion  for  the  electron  beam  are 
combined  together  with  the  partial  time  derivative  of  Eq.  (4.26)  to 
obtain  the  result 


+  2  +  + 
—  ij  -  o)  . ‘"Dp  +  u-  2  —  +  u-  —  pT 
bt  ^>2  2  0  0^^  bg 


-f  \ 


dz  \ 


(4.35) 


which,  for  convenience.  Is  written  In  the  following  form: 


1  -  1- 
dt  \  '  dt  ^2 


(I) 
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where 


[at 


(I) 


.at  '"o  az 


1  +  + 

2  S/b, 


fob  .+  .+ 

2  \az\ 


and 


- 

\at  ^2, 


2  ^  ^ 

D„  +  2  —  +  — 


(II) 

It  is  shown  in  Appendix  D  that 

2  2  r 


pb  2  “0 


at  ^  az 


2  2  1 

=  vii 


(I) 


2  1  i 


Ob 


hz 


and 


aS*  ao'  an*  'i  s  /  + 


hz  hz  hz 


bz 


Sd 

hz 


-\1 
J 


^at  S. 


(II) 


p'’'^  ■  (at^“°az)'’2 


2  _+ 

<0  +  u„€ 


f 

-  I  * 


O) 

-2£ 

0)  - 


1 

)  ] 

1  L 

aof  an"  a^D"  aS|' 


az  az 


az" 


az" 


r,  ^  /  aD"  an” 

=  ^  \  ^  az  ^  az 


We  know  also  from  Eq.  (4.27)  that 


1.-  1  -  i-fi  -  - 

at  P2  ‘  ""op  at  \  %t  \2  Vj>a 
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,(4.36) 


(4.37) 


(4.38) 


where  from  Eq*  (V.3I)  ve  have 


5  o  .  0)  ^ 

Pot,  —  V  =  ■*■  — P — 

^Op 

Together  Eqs.  (4.1^)  and  (4.11)  yield  the  expression 


(i^.39) 


By  substituting  Eqs.  (4.39)  and  (4.4o)  into  (4.38)  we  find  that 


(4.40) 


a  *  0) 

,  +  2  -_+ 

—  i  *0)  +  — 


a  P2  PP  2  2p^  0.2  I  1  1 


D:;  —  +  2D.  —  D 

^  az  ^  ^  az  ^ 


(i^.4i) 


From  the  condition  that  the  total  second  harmonic  current  (convection 
and  displacement  current)  is  zero,  it  follows  that 


:?r 


(lt.42) 


Substituting  in  the  above  expression  the  expressions  derived  before  for 
(a/at)i^  [Eqs.  (4.36)  and  (4.37)]  and  for  (a/at)!^  (Eq.  (4.4l)],  we 
obtain  tie  following  Inhomogeneous  differential  equation  for  the  second 


hairmonic  electric  field : 


^  t  -i  —  D 

Uq  bthz 


'2  ■  2  7^  3^^i  r 

\^2/ 


+  - ^  D^  —  .  D  —  D 


0)  ,  /  €-  \  d 

-i  Ip  •  D"  —  D" 

2Vo\^2/ 


(4.43) 


/  a^D:!'  aD’  aoil’N  cue  ^  /+^%' 

+  j-^(i  +  j)  D-  -4  -  U  j-^(|)  d;  ^  ^  , 

i«6^  ^  \  ^  8z2  az  az  /  i^€„  ^  \  ^  8z^  az  az  j 
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where  the  parameter  \i  ic  defined  to  be 


A 


2  2 
O)  CD  , 

,..PE. 


CD 


(4,44) 


For  the  case  in  which  there  is  no  plasma  (od  =  0  •  e-  =  1  ,  u  =  O)  the 
inhomogeneous  term  is  simplified,  and  we  have  instead  of  Eq,  (4.43)  the 
following  differential  equation: 


2 

Uq  '3t.^z 


(4  - 


£L 


2i  u 
0  0 


CD 


^0^0 


,  (4.45) 


12 

which  is  the  differential  equation  obtained  by  Paschke, 

Tlie  nonlinear  equation  [Eq.  (4.43)]  can  now  be  solved  for  the  second 
harmonic  field  by  substituting  the  solutions  for  the  first-order  terms 
into  its  right-hand  side.  To  solve  this  inhomogeneous  differential 
equation,  we  write  first  the  solution  of  the  associated  homogeneous 
equation: 


/  jh  z  -Jh  z 

.[A3.  ^  .B3. 


2j(ajt.p^z) 

e 


/  *  -Jhz  *  ^^2^ 

+  lA2e  +B2e 


-2j(aJt-P^z) 

e 


The  solution  of  Eq.  (4.43)  is  then  given  by  the  above  complementary 
function  added  to  a  particular  integral.  Appendix  E  contains  the  deri¬ 
vation  of  this  particular  solution  to  Eq.  (4.43). 
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With  the  condition  that,  at  z  =  0  ,  there  is  no  second  harmonic 
component  of  the  electric  field,  the  following  expression  is  obtained 
for  D*{z,t)  : 

^  i<o/oM\2  r  3(1  +  %) 

Dt(z,t)  =  —  *  <  (1  +  n) - 5 - ^  *  cos  2h.z 

HAv  L 


4(1  -  5^)  +  n(3  -  41^) 

+  - 5 -  cos  h  z  *  sin  2(cot-p  z) 

4r  -  1  ^  ® 


where 


sin  hgZ  *  cos  2(ajt-P^z 


(4.46) 


This  equation  can  be  conveniently  expressed  in  terms  of  Instead  of 

the  dimensionless  parameter  {  ,  with  the  help  of  the  following  relations: 
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where  the  term  proportional  to  neglected,  under  the  assumption 

that  h,/p  <  1  . 

1'  e 

The  definitions  of  and  have  been  given  before  and  are  re¬ 

peated  here  for  convenience; 

2  2 
0)  0) 

€  =  1  -  ^  ,  Gp  =  1  - 

^0)2  2  (2ao)^ 

These  expressions  can  be  considered  as  expressions  for  the  dielectric 
constants  exhibited  by  the  plasma,  at  excitation  frequencies  o)  and  2ud  , 
respectively.  Indeed,  it  can  be  shown  that,  for  a  infinite  plasma. 
Maxwell's  equations  can  be  written  in  a  simileur  form  to  that  applicable 
in  free  space,  but  with  a  dielectric  constant  given  by  or  ,  if 

we  assvune,  respectively,  time -dependence  for  the  rf  queuitities  of  the 
form  or  . 

The  solution  derived  above  for  the  second-harmonic  fields  presents 
some  interesting  features  which  will  be  the  subject  of  discussion  below. 

By  comparing  this  equation  with  the  corresponding  solution  derived  by 
Paschke,  in  the  absence  of  plasma,  one  can  see  that  the  presence  of 
plasma  introduces  a  significant  change  in  the  harmonic  amplitudes. 

The  parameters  u  and  i  which  appear  in  Eq.  (U.U6)  and  which  have 
been  defined  previously  as  being  given  respectively  by 


5  = 


> 


indicate  the  presence  of  plasma.  The  particular  values  u  *  0  ,  6  «  1 

express  the  absence  of  plasma  in  the  drifting  region  of  the  infinite  beam. 

In  this  case  Paschke*s  results  apply.  When  the  electron  beam  is  modulated 

at  exactly  the  plasma  frequency  (i.e.,  o)  =  oj  )  ,  then  we  find 

PP 


ji  =  00  ,  5  “  00  . 

The  resonance  at  O)  /2  is  attained  when  5*0 
PP' 


and 


)  • 
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Suppose  nov  that  the  beam  modulation  frequency  o)  is  gradually 

changed  from  the  condition  a)>a)  to  (x>  <  ca  /2  ,  We  shall  discuss 

PP  PP* 

the  Implications  in  the  second  harmonic  equation,  when  we  consider  each 
one  of  the  following  conditions: 


(i)  0)  >  05 

PP 


(li) 


05  >  05  >  05  /2 

PP  PP' 


(ill)  >  05 


Consider  first  (l);  05  >  05^^  .  In  order  to  compare  our  expression  for 

the  second  harmonic  field  with  that  of  Paschke  we  write  below  the 
corresponding  equation  derived  by  Paschke  for  a  velocity-modulated  beam 
of  finite  size,  with  no  plasma: 


05 


roM' 


-  i 


°  “pb  *  ^  ^ 


1 - 5 -  cos  2P  p(a5)z 

45^  .  1 


m  -  4^) 

+ - 5 -  cos  s  p(au5)z 

kc  - 1 


sin  2(o5t-P^z) 


(1^.50) 


where,  in  this  case,  the  dimensionless  parameter  i  denotes  the  ratio 
between  the  plasma  frequency  reduction  factors  for  the  finite  beam  at 
the  fundamental  and  second-harmonic  components,  respectively.  That  is, 
Paschke 's  parameter  S  is  5  *  p(«>)/p(2tt5)  .  If  the  following  corres¬ 
pondences  ere  then  established. 


p(o5)  -*■ 


1 

p(au)  -*■  -pr  , 

f2 

we  should  have  our  case  of  an  infinite  beam  in  eui  infinite  plasma,  the 
parameter  |  being  given  by  been  defined  before.  It 

is  evident  that  in  the  one-dimensional,  no  plasma  case  (05^^  *  O)  ,  we 
have  5  =  1  ,  as  stated  previously. 
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Consider  now  case  (ii);  <x>  >a)>(jD  /2,  When  the  above  condition 

is  satisfied,  it  may  be  seen,  from  Eqs.  (4.48)  and  (4,49)  that  becomes 
negative  and,  consequently,  the  argument  of  the  cosine  function  involving 

becomes  imaginary,  A  dependence  on  z  by  means  of  a  hyperbolic 
cosine  will  result,  and  this  means  that  a  growth  with  distance  occurs 
for  the  second  harmonic  fields. 

The  growth  of  the  fundamental  amplitude  which  exists  for  the  condition 
0)  <  ,  becomes,  in  this  case,  responsible  for  excitation  of  growing 

second  harmonic  fields.  If  the  beam  modulation  frequency  is  reduced  to 
a  value  such  that  ^pp/^  >  ^  y  only  is  a  negative  quantity  but 

6^  also  becomes  negative. 

It  follows  that  both  cosine  functions  in  Eq.  (4.48)  become  hyperbolic 
functions,  thus  indicating  growth  with  distance.  These  increasing  ampli¬ 
tudes  of  the  second  harmonic  waves  derive,  in  this  case,  from  nonlinearity 

due  to  the  growing  fundamented  waves  (o)  <  o)  )  and  also  from  interaction 

PP 

with  plasma  at  a  subharmonic  of  the  plasma  frequency  (a)^^2)  , 
a.  Second  Hetrmonic  Currents 

To  calculate  the  second  harmonic  plasma  current  density  we  use 
Eq.  (4.41)  to  obtain 


i*  =  (Cp  -  1) - I  •  (  D'  —  D  +  D  —  D" 

^2  2  at  ®  \  1  dz  1  1  dz  \ 

which  with  the  help  of  Eqs.  (4.23)  becomes 


,  (^.51) 


C  -  ('2  - 


^t 


\  0) 


o  /^l\ 

2  sin“  h^z  sin  2(a)t-P  z)  +  j  —  I  sin  2h^z  •  cos  2(a)t-P  z) 

®  \6  /  ® 

\  e/  J 


.  ('*.52) 
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The  second  harmonic  electron  beam  current  density  Is  derived  from 

the  equation,  1^  +1^  =  -  BDp/St  ,  where  use  Is  made  of  the  expression 

Pp  ^ 

for  the  plasma  current  density  [Eq.  (If. 52)].  Ety  neglecting  the  term 
proportional  to  have 


.  31^  \  ^  ^ 

1  *  -  €»  — -  +  - |"^1  1 —  I  ’  sln^  h,z  •  sin  2(ajt-p  z)  ,(^.53) 

^2  2  ^  2h)€^\a)/  \2  y  ^  ® 

or,  In  a  more  explicit  form,  with  the  help  of  Eq.  (4.46)  we  find 


(^-  foctt)  ^ 

‘“pb 


2n  23 

1  +  u  -  €^(l  +  -^)  COS  ^  Z 


-  ^1  +  (i  -  €j^(l  +  ^  2(<Dt-3^z)  ,  (^*5^) 

^2  J 

from  which  an  expression  for  ^  /2  could  be  easily  written, 

b.  Second  Hamonlc  Velocities 

The  electron  bean  second  harmonic  velocity,  as  derived  from 

Is  given  by  the  following  expression  (see  Appendix  F  for  details  In  the 
derivation) : 


o)/^  /cxm\^ 

-  -  2u^  -L-^  ~  • 


f 


6S(1  t  1 

l.i=  -  1  ■  5 


sin  2h^z  •  sin  2(a)t-3^z) 


4(1  -  5^)  +  ^(3  .  4^2 j 
45^  -  1 


sin  hgZ  •  sin  2(<ut-P^z) 


*2^  ^1*  *  2(<ot-3gZ)  >  , 
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where  terms  proportional  to  (h./P  )  and  (h^/P  )  have  been  omitted. 
This  is  consistent  with  the  assumption  that  second -harmonic  velocities 
are  not  excited  at  z  =  0  ,  as  it  may  be  seen  in  the  complete  expression 
for  v!^  ,  in  Appendix  F. 
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CHAPTER  V 
EXPERIMENTS 

5.1  DESCRIPTION  OF  THE  EXPERIMENTS 
a.  The  Beam-Plasma  Tube 

An  electron  beam-plasma  experiment  in  zero  magnetic  field  was  per¬ 
formed  with  a  tube  which  had  some  features  similar  to  the  tubes  described 
14 

by  Boyd.  The  essential  difference,  however,  lies  in  the  fact  that  we 
used  grids  for  modulation  of  the  electron  beam  in  an  attempt  to  investigate 
the  volume  space-charge  waves  and  to  study  the  behavior  of  the  beam  in  the 
plasma  under  saturation  conditions.  A  schematic  of  the  experimental  tube 
is  shown  in  Fig.  5>  and  a  photograph  of  it  is  also  included  (Fig.  6). 

Ibe  electron  beam  gun  employs  a  0.120"  diameter  planar  oxide  coated  cathode 
button  frcrn  a  2K28  Raytheon  klystron.  This  cathode  is  surrounded  by  a  heat 
shield  that  supports  the  accelerating  grid.  The  beam  is  velocity  modulated 
by  grids  which  were  taken  from  these  same  grldded  klystrons.  Identical 
grids  were  used  in  the  output  section  of  the  tube.  The  distance  between 
the  input  and  output  grids  is  approximately  10  cm.  This  region  consists 
of  a  12  mm  o.d.  nonex  glass  tube.  Beyond  the  output  grids  there  is  a 
beam  collector. 

A  helix  was  embedded  in  this  section  of  the  glass  tube  (0.383"  ) 

by  winding  0.010"  diameter  tungsten  wire.  The  helix,  which  is  nonsynchronous 
over  the  range  of  operating  beam  velocities,  acts  as  an  electrostatic 
shield  for  the  beam,  thus  avoiding  glass  charging  difficulties  and  per¬ 
mitting  cavity  measurements  of  plasma  density.  E:<periments  were  made 
earlier  v;lth  a  similar  tube,  one  idth  this  section  of  glass  coated  on 
the  inside  with  a  conducting  surface  (stannous  oxide). Although  gain 
was  observed  with  this  coated  tube,  measurements  of  plasma  density  by 
cavity  perturbation  techniques  were  found  to  be  unreliable  due  to  the  high 
screening  effect  exhibited  by  the  coating  surface.  No  magnetic  field  was 
used. 


'  ^It  was  because  of  this  geometry  that  the  calculations  in  Chapter 
III  assumed  a  metallic  boundary  outside  the  plasma.  It  is  not  obvious 
that  this  is  correct  for  the  stannous  oxide  film.  However,  as  is  shown 
in  Appendix  A,  the  nature  of  the  waves  is  not  markedly  altered  by  the 
assumption  of  a  dielectric  boundary;  for  the  nonsolenoldal  waves  not  at  all, 
and  for  the  solenoldal  v^aves  only  in  a  change  in  plasma  reduction  factors. 
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-Schematic  of  the  experimental  tiibe 


Pig.  6--Photograph  of  the  tube. 
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b.  Experimental  Procedure 

The  operating  frequency  was  selected  to  be  in  the  3000  Mc/s  region, 
usually  2890  Mc/s.  The  pair  of  cavities  used  in  the  tube  were  split 
cavities  designed  for  operation  in  the  10  cm  range  with  the  commercial 
gridded  klystrons,  from  which  the  grids  used  were  taken.  Tuning  was 
accomplished  by  tuning  screws  along  the  periphery  of  this  reentrant  cavity. 

To  obtain  quantitative  results  on  gain,  saturation  current,  etc., 
it  was  necessary  to  know  the  cavity  impedances.  From  a  perturbation  . 
measurement  an  accurate  value  of  R/Q  was  obtained.  The  result  was  R/Q 
=  57  ohms.  A  conventional  measurement  of  Q  ,  by  plotting  a  Q-clrcle, 
yielded  the  result  Q  =  130.  The  shunt  Impedance  of  the  cavity  was  then 
determined  yielding  R  =  Q(r/Q)  =  7^00  ohms. 

Bie  driving  source  was  a  1  milliwatt  1  kc/s  square-wave  modulated 
signal  generator.  When  it  became  necessary  to  drive  the  input  cavity 
at  large  signals,  this  signal  generator  was  fed  into  a  traveling  wave 
tube  amplifier  capable  of  delivering  1  watt  meiximum  output. 

Incident  power  to  the  input  cavity  was  measured,  using  a  directional 
coupler  and  rf  wattmeter,  llie  actual  power  absorbed  by  the  buncher  was 
obtained  from  the  knowledge  of  the  reflection  coefficient  in  the  Input 
line  as  determined  by  a  standing  wave  detector.  Then  from  the  knowledge 
of  the  rf  power  input  and  the  cavity  shunt  Impedance,  the  rf  voltage  at 
the  buncher  could  be  computed. 

Similarly,  from  the  knowledge  of  the  loaded  and  unloaded  Q  of 
the  output  cavity,  the  total  power  supplied  to  the  output  gap  could  be 
calculated  in  terms  of  the  measured  power  reaching  the  external  load. 

The  beam  coupling  coefficient  M  Is  given  in  terms  of  the  dc  gap 
transit  angle  d  »  where  -t  is  the  gap  width.  This  distance  -t 

between  the  grids  constituting  the  gap  was  measured,  yielding  the  result 
-t  -  O.okO".  It  follows  then  that  M  =  sin(d/2)/(<i/2)  value  M  =  0.9» 

for  a  beam  velocity  corresponding  to  500  v  and  frequency  of  2890  Mc/s  . 

The  output  cavity  was  heavily  overcoupled  and  had  a  coupling  coefficient 
P  s  Q^Qext  equal  to  5-  The  rf  beam  cvirrent  at  the  output  cavity  was 
then  calculated  in  terms  of  the  known  measured  quantities  by  using  the 
expression 

^  I  \ ^  2(1  + 

(ij  ■ 
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where  denotes  the  power  to  the  external  load,  and  is  the  dc 

beam  current.  It  v/as  found  unnecessary  to  make  measurements  for  different 
values  of  the  beam  current,  since  the  effect  of  beam  loading  on  the  output 
cavity  was  negligible.  In  the  final  part  of  the  experimental  work,  to  be 
described  later,  attention  was  concentrated  upon  the  possibility  of 
harmonic  generation.  A  nev;  input  cavity  v/as  designed  and  built  to 
operate  in  L-band,  and  to  be  used  as  the  buncher.  The  fundamental 
operating  frequency  was  chosen  to  be  1770  Mc/s  and  the  output  cavity 
was  tuned  to  twice  this  frequency. 

c.  Electron  Beam 

(1) .  Focusing  and  Interception  Problems 

It  v/as  found  difficult  to  achieve  perfect  focusing  of  the  electron 
beam  along  the  axis  of  the  tube.  It  was  in  trying  to  fulfill  this 
idealized  requirement,  with  the  arc  off  and  no  magnetic  field,  that  trouble 
was  first  encountered. 

Nevertheless,  after  proper  voltage  adjustments,  it  was  possible, 
with  no  plasma^  to  determine  that  the  beam  v/as  reasonably  focusing  in  the 
helix  portion  of  the  tube  by  visual  observation  of  the  luminosity  it 
produced  along  its  path. 

The  dc  behavior  of  the  electron  beam  in  this  region  is  controlled 
by  the  accelerating  voltage  impressed  on  the  input  gap  and  the  dc  voltages 
applied  to  the  helix  and  output  grids,  respectively.  With  the  plaama  on, 
the  beam  could  not  be  seen.  We  re]y,in  this  case,  on  the  positive  ion 
neutralization  due  to  the  plasma  ions  which  should  improve  the  beam 
focus  somewhat.  The  dc  beam  current  reaching  the  output  region  can 

be  considered  as  a  measure  of  the  degree  of  this  improvement. 

The  beam  transmission  to  the  collector  was  measured  to  be  of  the 
order  of  30^*  Interception  by  the  focusing  grid  near  the  cathode  and 
the  two  pair  of  grids  v/hich  form  the  gaps,  besides  collisions  v/ith  neutral 
molecules,  are  responsible  for  this  relatively  poor  transmission. 

(2) .  Electron  Beam  Parameters 

For  the  reasons  discussed  above,  an  accurate  measurement  of  the 
diameter  of  the  electron  beam  is  not  possible.  We  assumed  in  the  cal¬ 
culations  made  that  the  beam  diameter  in  the  interaction  path  is  approximately 
the  diameter  of  the  cathode  button. 
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The  beam  current  Is  partially  controlled  by  applying  to  the  focussing 
grid  positive  (or  negative)  voltages  with  respect  to  the  cathode.  Un¬ 
fortunately,  it  was  found  difficult  to  maintain  a  constant  beam  current, 
as  measured  at  the  cathode,  over  a  wide  range  of  beam  voltages.  Arcing 
between  the  gun  electrodes,  due  to  ions  in  the  gun  region,  often  occurred, 
therefore  limiting  the  operation  of  the  beam  gun.  Moreover,  the  beam 
cxirrent  was  shown  to  vary  as  much  as  20^,  as  the  arc  current  was  varied 
in  the  course  of  the  experiment.  Positive  ion  bombardment  of  the  emitting 
surface  covild  not  be  avoided. 

By  recording  the  magnitudes  of  the  dc  collector  current  for  different 
plasma  currents,  and  by  comparison  with  the  readings  at  the  collector 
when  only  the  plasma  is  on,  the  actual  beam  collector  current  is  recorded. 
In' the  present  experiments,  the  electron  beam  voltage  .anges  from  100 
to  500  volts.  Values  of  beam  cathode  ciirrents  of  1  ma  and  1.5  ma  were 
convenient  for  our  purposes.  The  beam  plasma  wave  number  for  a 

beam  voltage  500  v,  and  a  beam  current  1  ma,  3  nm  diameter  is  such  that: 
Pp^a  =  0.05,  where  a  is  the  beam  radius.  At  a  frequency  of  2890  Mc/s 
and  a  beam  velocity  corresponding  to  500  volts,  the  beam  electronic  wave 
number  p  --  ciVu.  is  such  that  p  a  -  2  . 
d.  Plasma 

(1).  Characteristics  of  the  Arc  Discharge 

The  plasma  used  in  these  experiments  consists  of  the  positive 
column  of  a  low  pressure  hot  cathode  merevury  arc  discharge.  The 
theoretical  calculations  of  beam  interaction  in  earlier  chapters  assvimed 
a  completely  quiescent  ("cold")  plasma  with  no  collisions.  To  estimate 
the  effects  of  departiures  from  this  ideal  condition  in  later  sections 
we  shall  need  some  values  of  the  plasma  temperature  and  colllson  fre¬ 
quency. 

The  vapor  pressure  of  mercury  at  room  temperature  (26.8  C)  is 
2.1  X  10~^  mm  Hg.  The  equivalent  plasma  electron  traperature  (T^)  is 
defined  in  terms  of  the  random  energy  of  electrons  according  to 


-  aw  ^  =  1  kT 
2  t  2  e 


(5.1) 


where  v^  is  the  average  thermal  velocity  of  the  plasma  electrons  and 
k  is  the  Boltzmann  constant. 
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Boyd  reported  measurement  of  this  plasma  electron  temperature 
with  a  Langmuir  probe.  An  average  value  of  =  36OOO  K  was  obtained 
in  a  tube  of  internal  diameter  equal  to  1.0^  cm.  The  straight  part  of 
our  tube  where  the  beam-plasma  interaction  takes  place  has  an  inside 
diameter  of  0.9  cm.  We  may  then  assume  this  value  of  T^  as  approximately 
correct  for  our  case.  A  random  energy  for  the  plasma  electrons  of 
approximately  U.7  ev  woiild  correspond  to  this  plasma  temperature. 

The  collision frequency  as  calculated  from  the  mean  free  path  for 

lU  , 

collisions  of  electrons  with  neutrals  was  approximately  37  Mc/s. 

Collision  effects  are  therefore  negligible  for  this  value  of  the  collision 

frequency  at  the  modulation  frequencies  of  3000  Mc/s  or  I800  Mc/s  used 

in  the  experiments. 


(2).  Plasma  Electron  Density;  Cavity  Perturbation  Measurements 

The  plasma  electron  density  is  easily  calculated  from  the  formula 


f 


2rt  V 


(5.2) 


where  n  designates  the  number  of  electrons  /csur  ,  where  f^  is  the 

plasma  frequency  of  the  plasma,  is  the  permittivity  of  free  space, 

e  and  m  the  charge  and  mass  of  the  electron,  respectively.  At  a 

plasma  frequency  of  3000  Mc/s,  the  plasma  electron  density  turns  out  to 
XX  3 

be  1.12  X  10  electrons/cm^  .  It  is  convenient  for  this  purpose  to 

15 

consider  the  plasma  as  a  charge-free  region  exhibiting  a  frequency- 

2  2 

dependent  dielectric  constant  given  by  c/c^  =  1  -  0)^^  /o)  ,  (when  no 
magnetic  field  is  present  )  and  one  can  measure  the  dielectric  constant 
by  perturbing  the  resonant  frequency  of  a  cavity  by  introducing  a  plasma. 

A  split  cylindrical  cavity  supporting  a  mode  was  designed 

for  this  purpose.  Measurements  were  made  of  the  change  of  resonant 
frequency  of  this  mode  when  the  plasma  density  was  varied  by  means  of 
the  arc  current.  The  frequency  deviations  allow  us  to  calcvilate  the 
average  electron  plasma  density  with  the  help  of  a  perturbation  formula 
Figure  7  is  a  plot  of  the  square  of  plasma  frequency  against  arc  current. 
Indicating  approximately  the  expected  linear  dependence  of  the  plotted 
variables . 
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5.2  EXPERIMENTAL  RESULTS;  DISCUSSION 

Power  measurements  were  taken  both  at  the  input  and  output  of  the 
tube,  according  to  the  procedure  outlined  in  the  previous  section.  First, 
the  electron  beam  was  modulated  vdth  a  low  power  signal  and  measurements 
were  taken  of  the  output  power,  for  a  fixed  input,  as  a  function  of  the 
plasma  density.  The  plasma  density  was  changed  by  varying  the  arc 
current.  Figure  8  is  a  plot  of  the  output  power  versus  arc  current  with 
a  low-level  input  signal  (  input  power  1  mwatt).  By  comparing  this  curve 
with  the  data  showing  the  plasma  frequency  dependence  on  arc  current 
(Fig.  7)>  one  can  see  that  at  the  operating  frequency  of  2890  Mc/s,  the 
maximum  output  signal  appears  near  the  expected  value  for  resonance 
(i.e.,  f  ~  fp).  Although  no  net  gain  was  observed,  it  may  be  seen  that 
this  maximum  output  corresponds  to  a  gain  due  to  the  plasma  of  approximately 
5.5  db.  With  the  same  dc  conditions  for  the  electron  beam  (beam  velocity 
and  beam  current  unchanged),  this  experiment  was  repeated  with  a  larger 
driving  power  (20  mwatts)  which  was  the  order  of  power  necessary  to  cause 
rf  saturation  of  the  beam.  A  plot  of  the  output  power  as  a  function  of 
arc  current  is  given  in  Fig.  9-  It  is  seen  that  the  saturation  power  is 
about  10  db  above  the  low-level  power,  for  an  increase  in  drive  of  13  <ib, 
and  so  the  saturation  "gain''  is  3  db  below  the  small-signal  gain.  Similar 
data  taken  with  the  same  beam  velocity  and  modulation  frequency,  and  20 
mlllii;atts  rf  input  power,  but  with  a  beam  current  of  1.5  ma  instead 
of  1  ma,  are  plotted  in  Fig.  10.  For  this  value  of  input  power,  the 
gain  due  to  the  plasma  is  7*5  db  at  an  arc  current  of  approximately  220 
ma,  which  corresponds  to  an  average  plasma  density  such  that  the  resonance 
condition  is  nearly  satisfied.  It  is  interesting  to  notice  in  these 
curves  that,  as  we  might  expect  from  the  theory,  for  values  of  the  arc 
cvirrent  above  that  for  maximvun  gain,  the  output  signal  falls  less  rapidly 
as  current  is  changed  than  it  does  for  arc  currents  below  this  value. 

The  broad  shape  of  the  curves  might  be  explained  as  a  consequence  of 
inhomogeneity  in  plasma  density  and  as  a  result  of  the  poor  focusing 
of  the  beam.  As  pointed  out  before,  the  electron  beam  was  not  as  well 
focused  along  the  axis  of  the  plasma  column  as  would  be  desirable. 
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-Output  power  versus  arc  current 
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Fig*  9 — Output  power  versus  arc  current  for  input  saturation  power 


Fig.  10 — Output  power  versus  input  power  for  different  values  of  arc 
currents  (beam  current  ~  l.p  raa). 
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The  interaction  process  may  tcke  place  there  at  different  distances  from 
the  axis  of  the  plasrar,  which  is  not  uniform  over  its  cross  section. 
Nonuniformity  in  the  plasma,  i.e.,  only  its  length,  was  presumably  one 
of  the  main  effects  contributing  to  the  small  values  of  gain  which 
were  measured.  Actually,  it  wcs  observed  that  the  plasma,  did  not  completely 
fill  the  glass  tube  along  the  pc?th  of  the  electron  beam,  from  input  to 
output  cavity.  Regions  near  the  cavities  were  found  to  be  free  of  plasma. 
Under  these  circumstances,  the  electron  beam  leaves  the  plasma  region 
and  drifts  in  a  plasma-free  space  with  a  length  of  approximately  2  cm, 
before  entering  the  output  cavity.  This  corresponds  to  a  length  =  *7^ 

and  such  a  drift  region  may  reduce  the  gain  and  maximum  rf  current. 

An  appreciable  reduction  in  the  output  levels  would  also  be  expected 
by  virtue  of  thermal  effects  and  collisions  in  the  plasma,  which  we  shall 
discuss  below.  However,  the  calculations  will  show  that  at  maximum  these 
effects  cannot  completely  account  for  the  low  gain  and  it  is  probable 
that  the  nonuniformity  of  the  plasma  is  the  principal  source  of  gain 
reduction. 

There  are  several  interesting  aspects  to  the  measurements  of  the 
output  power  versus  input  power  v/hich  are  plotted  in  Figs.  11  and  12. 

These  measurements  were  taken  successively  at  different  arc  currents 
whose  values  correspond,  respectively,  to  the  conditions  at  which  the 
modulation  frequency  is 

(a)  below  the  plasma  frequency 

(b)  near  the  plasma  frequency 

(c)  above  the  plasma  frequency. 

Power  measurements  with  only  the  electron  beam,  in  the  cbsence  of  plasma, 
are  slso  included.  It  is  important  to  note,  from  these  curves,  the 
significant  increase  in  the  saturation  output  power  when  the  plasma  is  on. 

Tlie  measurements  given  in  Figs.  11  and  12  are  replotted  in  Figs.  13 
and  1^  in  terms  of  the  output  current  versus  input  voltage.  In  these 
graphs,  the  dc  beam  current  given  is  that  baaed  on  collector  current  rather 
than  total  eraitted  current. 

It  is  known  that  the  simple  first-order  bunching  theory  due  to  Webster 
predicts  for  a  two-cavity  klystron  that 

=  2IqJ^(X)  ,  (5.3) 
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Fig.  11 — Output  power  versus  input  power  for  different  values  of  arc 
current  (beam  current  ~  1  ma). 
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Pig.  13 — Pundameutal  component  of  the  current  versus  the  bunching  parameter  (beam  current  ~  1  ma). 
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where 


=  rf  beam  current 
Iq  =  dc  beam  current 
X  =  nNQ  =  bunching  parameter  , 

derived  for  a  beam  of  infinite  extent^  without  any  consideration  of  space- 
charge  effects. 

From  the  observed  output  currents^  with  no  arc  present,  it  may  be  seen 
that  the  maximum  value  of  the  rf  current  is  approximately  0.8  ,  at  a 

value  of  the  parameter  approximately  1.9;  this  value  for  X  is  indeed  near 
the  theoretical  value,  1.84,  for  which  the  maximum  of  J^(X)  would  be  ex¬ 
pected  to  occur.  The  fact  that  the  observed  value  of  the  maximum  current 
is  below  l.l6,  as  it  should  be  in  the  simple  formula,  would  indicate  the 
effect  of  the  space-charge  debunching  of  the  electron  beam  which  drifts 
along  a  tube  of  length  L  such  that  ^  •  However,  it  should 

also  be  noted  that  the  observed  minimum  of  the  current  is  at  a  value  of  the 
bunching  parameter  X(X  =  3.6),  very  close  to  the  value  of  the  argument  of 
the  Bessel  function  which  corresponds  to  a  zero  of  this  function 

(X  =  3.83). 

The  discrepancy  observed  at  small  signal  levels  between  the  theoretical 
values  for  the  current  as  given  by  the  Bessel  function  dependence  and  the 
corresponding  experimental  values  may  also  be  accounted  for  by  the  finite- 
beam  space-charge  effects. 

When  the  argument  of  the  Bessel  function  is  small,  according  to  the 
Ramo  single-mode  approximation,  the  expression  for  current  in  a  finite 
beam  is  simply 


I 


1 


I 


0 


(5.^) 


In  the  above  equation  represents  the  "effective"  or  "reduced"  plasma 

frequency  of  the  beam. 

A  similar  expression  applies  to  the  case  of  a  beam  passing  through  a 
plasma  with  the  reduction  factor  being  due  to  the  presence  of  the  plasma. 

To  make  any  comparison  with  experiments  for  this  case  (beam-plasma),  we 
must  first  Justify  an  approximation  made  in  Chapter  III.  In  considering  the 
nonsolenoidal  waves  the  assumption  was  made  that  h/p^  «  1  .  This  obviously 
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in  not  true  at  6o  =  o)  vhorc  the  Gitnplo  theory  'j;ive.n  an  infinite  value 
for  this  ratio..  Tlic  arjsuraption  can  bo  Justified  even  for  this  value  of 
(0  if  one  considers  the  effects  of  thermal  velocities  of  the  plasma  electrons. 
The  effect  of  chese  thermal  velocities  has  been  calculated  for  a  plasma- 
beam  system  of  infinite  extent  but  the  calculations  are  applicable  to  the 
nonsolcnoidal  waves  of  a  finite  system,  since  the  difficulties  with  h/P^ 
arise  in  the  same  form  for  both  cases. 

For  the  electron  beam  and  the  plasma,  both  of  infinite  extent,  it 
has  been  shown  [Chapter  III,  Eq.  (3.127)]  that  the  fundaraental  coraponent 
of  the  current  is  Given  by  the  expression 


ii'V. 
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(5.5) 


D  '  j 

Comparing  the  above  two  expressions,  (p.^)  and  (>#5),  the  analogy  at  small 
.signals,  between  the  finite-beam  case  and  tlio  infinite  beara-plasma  case, 
is  apparent.  The  correspondence  is  made,  as  in  Chapter  IV,  between  the 


CO 


case  and  l/^ 


plasma  frequency  reduction  factor  p 
in  the  beam-plasma  case. 

Plo^/ever,  while  p  is  always  less  than  unit,’',  as  already  stated, 

the  factor  could  he  allov/ed  theoretically  to  assurae  values  greater 

han  unity,  to  boccrae  iraaginary  when  co  <  ,  and  even  to  be  inl’inite, 

at  CO  =  CO  ,  in  a  loccloss,  unh’or.n  placraa.  The  ideal  situation  of  a 
PP 

lossless  plasma  is  not  encountered  in  practice.  In  our  present  experiment, 
it  may  be  shown  that  with  the  inclusion  of  thermal  effects,^ for  the  t^^plca 
bcara  parameters  used  here,  this  factor  turns  out  to  be  greater  than  unity 
but  not  infinite  at  the  plasma  resonance  frequency. 

Tlie  dispersion  relation  for  the  infinite  bca:n  In  a  plasma  has  been 
col/ed  nuinerically  by  taking  into  account  the  tlicraial  raotion  of  t.:e  plasma 
^•lectrons  and  collisions  (see  ref.  l4).  At  co  =  co^^  it  was  found  that, 
neglecting  collisions, 


PP 


(5.6) 


TlI 


It  can  be  cho^m  that  collieion  effects  are  nrech  smaller. 
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where  R  denotes  the  ratio  of  the  beam  electron  energy  to  the  averace 
random  ener'jy  of  the  plasma  electrons,  l.e., 


R 


1  2 
2 

I  kT 


(5.7) 


Assuming  a  500  volt  electron  beam,  and 
For  a  beam  ciu^rent  of  1  ma,  we  have  at 
3  liMc/s: 


T  =  36000°  K,  we  have  R  ~  100  . 
e 

0  =  0)  ,  for  a  frequency  equal  to 


^  -4 

6  X  10 


It  foil  ws  then  that 


which  justify  the  approximations  made  <  i)  in  the  analysis  of  the 

problem  studied  in  Chapter  III.  The  nonsolenoidal  solutions  were  shown 
in  that  chapter  to  yield  the  following  expression  for  the  volume  beam 
current : 


I 

v 


2  h. 

j. 


sin  (hs)  +  •  •  • 


) 


applicable  for 
a)  approaches 


Cl)  <  ,  which  should  be  compared  with  Eq.  (5*  5)*  When 

^PP  ^  using  the  quantities  calculated  above,  v/e  have 


^  10  ,  h  1.3  nepers/cm  . 

Under  these  conditions,  the  growth  constant  for  these  nonsolenoidal  cxirrcnts 
is 

g  =  0.68  h.  =  8.63  X  1.3  =  11.2  db/cm 
n.s.  1  ' 

Tlic  solenoidal  currents  which  are  associated  with  pura.l;;,"  surface  wave  prop- 
a:;ation  have  also  been  investigated  in  Chapter  III,  since  they  also  constitute 
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peu:t  of  the  excitation  produced  by  grid  coupling  to  the  electron  beam. 
However,  as  was  shown  before,  these  waves  would  lead  to  rates  of  growth 
which,  as  compared  to  the  nonsolenoidal  waves,  are  reduced  by  a  factor 
F  due  to  the  finiteness  of  the  system.  This  factor  was  calculated  as  a 
function  of  the  ratio  of  the  plasma  and  electron  beam  diameters,  and  also 
of  P^a  ,  with  the  assumption  that  ^  •  This  reduction  factor  is 

obtained  from  the  curves  plotted  in  Fig.  2.  In  our  case,  we  shall  assume 
that  we  may  take  account  of  thermal  velocity  in  the  same  way  as  for  the 
infinite  plasma  and  still  use  the  same  reduction  F  to  estimate  the  gain 
of  the  finite  system.  Therefore,  by  assuming  b/a  =  3  ,  =  2  ,  we  have 

F  s  0.6  . 

It  follows  therefore  that  the  growth  constant,  for  these  solenoidal  currents, 
at  the  plasma  resonance  frequency,  should  be 

g  =  Fh^  =  0.6  X  1.3  =  0.78  nepers/cm  , 
or 

g  =6.9  db/cm  . 

®1 

These  theoretical  results  indicate  that  there  are  no  difficulties 
due  to  infinities  in  the  gain  chsuracteristlc.  It  should  be  pointed  out, 
however,  that  the  small  signal  gain  calculated  here  for  both  kinds  of  waves 
is  much  larger  than  actually  measured.  As  previously  stated  this  discrep¬ 
ancy  is  probably  due  to  inhomogeneities  in  the  plasma  and,  in  peurticuleu*,  the 
existence  of  a  short  plasma-free  region  Just  preceding  the  output  cavity 
could  affect  the  small  signal  gain  greatly.  Additional  evidence  for  this 
is  shown  by  the  small  decrease  in  gain  for  Increase  in  plasma  frequency 
above  resonance  (400  ma,  ac  current,  see  Figs.  11,  12).  For  a  homogeneous 
plasma  density  the  difference  between  250  ma  euid  400  ma  should  be  much 
greater . 

a.  Satiu*ation  Effects. 

Under  large  signal  saturated  conditions,  with  crossing  trajectory, 
a  space-cheurge  wave  theory  is  never  valid,  and .  one  can  only  consider  the 
experimental  results. 
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An  investigation  of  the  experimental  curves  shows  that,  neeu:  resonance, 
the  rf  output  currents  reached  maximvmi  values  around  1.3  Iq  and  1.5  Iq  , 
which  respectively  correspond  to  total  beam  cathode  currents  of  1  ma  and 
1.5  raa.  The  curves  were  calculated  on  the  basis  of  the  dc  beam  current 
reaching  the  collector  which  was  found  to  be  about  30  ^  of  the  total  beam 
ciirrent.  These  values  for  the  rf  cxurrents  should  be  compared  with  the 
maximum  theoretical  value  for  the  electron  beam  in  a  drift  tube,  in  the 
absence  of  plasma.  As  was  shown  before,  in  theory  the  fundamental  component 
of  the  beam  current,  with  no  plasma,  heis  a  maximxam  for  ~  * 

The  reasons  for  the  Increase  in  the  ratio,  l/l^  ,  in  the  presence 
of  the  plasma  as  con^jared  to  this  ratio  in  the  absence  of  the  pleisma  can 
be  understood  by  considering  the  nature  of  the  interaction  mechanism  between 
the  electron  beam  and  the  plasma  at  frequencies  vrtiere  gain  occurs.  As  can 
be  seen  from  the  theory,  for  o)  <  the  plasma  exhibits  a  negative  di¬ 

electric  constant  and  in  its  response  to  the  current  in  the  electron  beam 

behaves  like  an  inductive  medium.  For  this  condition  O)  <  o)  the  modulated 

PP 

beam  is  continuously  peissing  through  an  inductive  mediimi.  It  is  known 
from  the  theory  of  multi-cavity  klystrons  that  if  the  cavities  close  to 
the  output  end  eu*e  inductively  detuned  then  the  resultant  modulation  pro¬ 
duced  by  these  cavities  is  such  as  to  distort  the  bunching  in  the  electron 
beam  in  a  way  so  as  to  produce  larger  rf  currents  than  if  one  had  merely  a 
beam  drifting  through  a  drift  tube. (or  through  cavities  tuned  to  resonance). 
This  is  known  theoretically  emd  experimentally,  and  it  has  been  well  estab¬ 
lished  that  detuning  (as  is  usually  the  case)  the  second  to  last  cavity  in 
a  klystron  will  Increase  the  efficiency  by  increasing  the  saturated  value 
of  the  rf  current,  beyond  idiat  is  possible  with  a  series  of  cavities  all 
tuned  to  resonance.  Presumably,  the  same  phenomena  occxirs  in  the  case 
described  here  in  the  plasma  region  near  the  output  cavity.  In  this  region 
the  rf  cvurrent  in  the  beam  produces  fields  in  the  plasma  such  eis  to  produce 
better  bunching  than  if  the* medium  were  either  not  present  or  not  Inductive. 
The  Increase  In  the  rf  current  raeasiured  here  Indicated  this  phenomena  Is 
taking  place. 

The  reasons  for  the  discrepancy  in  the  ratio  of  rf  to  dc  for  the  two 
values  of  dc  current^  shown  in  the  figures,  is  not  easily  understood  and  ma;y  oe 
due  to  experimental  inaccuracy,  e.g.,  dlscr^ancies  between  beam  cxirrent 
measured  at  the  collector  and  the  true  cvirrent  effective  in  the  Interaction. 
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b.  Second  Harmonic  Experiments > 

Experiments  similar  to  the  ones  described  above  for  measuring  second 
harmonic  currents  were  performed  by  using  a  different  pair  of  cavities  ap¬ 
propriate  for  detection  of  second-harmonic  signals.  The  input  cavity  was 
tuned  at  1770  Mc/s  and  the  output  cavity  at  twice  this  frequency.  Measure¬ 
ments  of  output  power  at  second  harmonic  versus  arc  current  are  plotted  in 
Fig.  15.  The  observed  peaks  at  values  of  arc  current  around  90  ma  and 
290  ma  are  found  to  correspond  to  the  conditions  that  the  input  frequency 
0)  is  o)  O)^  and  co  ^  >  respectively,  at  a  modulation  frequency  equal 

to  1770  Mc/s. 

Figure  16  shows  results  of  input  and  output  measurements,  indicating 
a  maximum  output  power  for  the  second  harmonics  at  00  ~ 
analysis  made  in  a  previous  chapter,  for  infinite  beam-plasma,  it  may  be 
seen  that  maximum  values  for  the  second  harmonic  output  should  be  expected 
for  o)  =  ODp  and  cu  =  ^0^/2  .  The  experiment  confirmed  the  theory  by  re¬ 
vealing  that,  for  a  fixed  input  frequency,  these  maxima  occur  at  a  sub¬ 
harmonic  of  the  plasma  frequency  (0)^/2)  and  also  at  the  plasma  frequency 

0)  . 

P 

The  experimental  data  were  also  plotted  in  a  normalized  fashion, 
similar  to  the  plots  of  Figs,  ih  and  15,  so  that  they  could  be  cofl^ared 
directly  with  the  results  obtained  previously  for  the  fundamental  components. 

The  plot  shown  in  Fig.  17  of  the  second  harmonic  ciirrent  at  the  out¬ 
put  cavity  with  respect  to  twice  the  bunching  parameter  X  ,  calculated  at 
2890  Mc/s,  was  derived  from  the  observed  curves  of  input  power  versus  out¬ 
put  power.  The  low  R/Q*s  of  the  input  and  output  cavities  used  limited 
the  power  levels  obtainable.  Nevertheless,  one  can  see  the  significant 
effect  due  to  the  plasma  when  comparison  is  made  with  the  situation  Mhere 
only  the  electron  beam  is  present.  The  saturated  output  current  is  shown 
to  increase  by  a  factor  of  approximately  two  when  the  arc  current  is  changed 
from  zero  (no  plasma)  to  ®  290  ma  (o)  o) p/2).  The  fundamental  com¬ 

ponents  were  not  measured  at  the  same  time  since  it  was  not  possible  to 
tune  the  output  cavity  to  a  frequency  as  low  as  1770  Mc/s.  However,  we 
may  make  some  comparison  between  the  relative  magnitudes  of  fundamental  com¬ 
ponents  of  current  and  second  harmonic  by  considering  some  of  the  experimental 
curves  that  were  taken  at  2890  Mc/s.  For  the  latter  in  the  saturation  region 
(input  power  of  30  mw,  for  instance),  the  maximum  value  of  a  fundamental 
component  of  current  was  approximately  1.3  Iq  ,  whereas  in  the  same  range 
of  frequency,  the  maximum  value  of  current  as  a  second  harmonic  component 
was  about  0.2  . 
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volta^^je  =  50J 


100 


Fig.  1, '"--Second  harmonic  component  of  the  current  versus  2X  =  2irK<^ 
(calculated  at  2690  Mc/s). 


5.3  SUMMARY  AI®  CONCLUSIONS 

The  experimental  data  described  in  this  chapter  agree  qualitatively 
with  the  predictions  of  the  theory.  We  found  that  the  gain  measured  be¬ 
tween  two  cavities  was  not  appreciably  affected  by  the  presence  of  a  plasma 
until  the  plasma  density  reached  a  value  approximately  equal  to  the  modulating 
frequency.  For  this  condition,  an  increase  in  gain  was  noticed  although 
numerically  this  gain  due  to  the  plasma  was  considerably  less  than  would 
be  indicated  by  the  theory. (even  with  thermal  velocities  in  the  plasma). 
However,  there  were  imperfections  in  the  plasma  column,  principally  in¬ 
homogeneities,  and  an  actual  short  region  preceding  the  output  cavity  which 
was  free  of  plasma,  all  of  which  would  tend  to  drastically  reduce  the  ex¬ 
pected  gain. 

The  interesting  feature  of  the  small  signal  performance  however  was 
that,  at  plasma  densities  corresponding  to  higher  plasma  frequencies  than 
the  modulating  frequency,  one  still  showed  gain  as  compared  to  a  no  plasma 
case.  In  agreement  with  the  qualitative  prediction  of  the  theory,  there 
was  a  marked  asymmetry  in  the  gain  shown  for  densities  lower  than  plasma 
frequency  as  compared  to  above  plasma  frequency. 

Possibly  the  most  interesting  feature  of  the  experimental  results  was 
the  saturation  current  measured  in  the  presence  of  a  plasma  with  the  plasma 
frequency  approximately  equal  to  the  modulati:**g  frequency.  Under  these 
circumstances  there  was  a  marked  increase  in  the  saturation  current,  as 
compared  to  the  no  plasma  or  low  density  plasma  cases,  l^pical  niombers 
for  the  ratio  of  rf  current  to  the  dc  current  ranged  from  1.3  to  1,5  in 
the  case  >^ere  the  plasma  frequency  was  approximately  equal  to  the  modulating 
frequency,  whereas  this  ratio  was  about  0.8  for  the  no  plasma  case.  This 
behavior  is  attributed  to  the  distorted  bunching  effects  of  an  inductive 
medium  acting  on  an  electron  beam  and  should  have  interesting  implications 
for  the  use  of  the  plasma  direction,  in  an  actual  operating  device.  Measure¬ 
ments  of  the  heurmonic  content  of  an  electron  beam  passing  through  a  plasma 
showed,  as  expected,  marked  effects  when  the  plasma  frequency  was  approximately 
equal  to  the  modulating  frequency  or  twice  the  modulating  frequency.  Biese 
results  may  have  implications  for  using  a  plasma  device  as  a  frequency 
multiplier. 
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APPENDIX  A 


ELECTRON  BEAM  IN  A  PLASMA;  DIELECTRIC  WALLS 


The  theoretical  model  which  was  considered  in  Chapter  III  does  not 
completely  apply  to  our  experimental  tube.  In  this  tube,  there  are  no 
conducting  v/alls  surrounding  the  plasma  colimin.  Dielectric  walls  (glass) 
should  be  assimied  Instead. 

In  this  appendix  a  brief  discussion  is  given  which  shows  that  the 
presence  of  metr.llic  walls  is  not  relevant  to  the  validity  of  the  con¬ 
clusions  stated  in  Chapter  III  concerning  the  fields  vanishing  in  the 
region  outside  the  beam.  Refer  to  Fig.  1,  and  assume  that  In  the 
dielectric  region  (region  III),  the  longitudinal  electric  fields  vary 
as 


E 

^III 


CH^^^^Pr)  . 


As  before,  in  region  II,  we  have 


The  boundary  conditions  at  r  =  a  are 


and 


E  (a)  *  E  (a)  (A.l) 

H  ^11 


r*a 


0 


(A.2) 


which  results  from  £q.  (3.^2),  by  assuming  nonsolenoldal  solutions  (propa¬ 
gation  constants  3  =  1  h)  Inside  the  beam  region. 

At  r  =  b  (plasma-dielectric),  we  have 


E  =  E  .  (A. 3) 

"II  "III 
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It  is  apparent  that  Ag  =  Bg  *  C  *  0  is  a  solution  of  the  system  of 
equations  represented  by  (A-l),  (A-2),  (A-3).  This  solution  is  com¬ 
patible  with  the  nonsolenoidal  fields  assumed  inside  the  beam.  It  is 
known,  from  the  analysis  in  Chapter  III  with  a  metallic  boundary,  that 
these  nonsolenoidal  solutions  in  region  I  are  determined  without  any 
effects  imposed  by  the  region  outside  the  beam.  The  treatment  of  this 
appendix  merely  shows  formally  that  the  same  holds  for  an  unbovnded 
dielectric  region  outside  the  plasma. 


APPENDIX  B 

CALCULATION  OF  THE  TOTAL  VOLUME  CURRENT 

To  calculate  the  volume  current  density,  ve  consider  the  separate 
contributions  of  the  electron  beam  and  plasma,  and  in  both  cases  ve 
take  into  account  the  nonsolenoidal  and  solenoidal  solutions  for  the 
fields. 

1.  Electron  Beam 

a.  Nonsolenoideil  component  of  the  beam  current  density. 

We  have: 


zb 


-Jhz  i[l'> 

e  +  - 


Jhz 


- LLi  cos(hz)  -  J  -25 - 2^)sln(hz)  . 


But,  we  know 


zb  zb 


b  '  ^  b 
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u. 


and 


l(l)  .  l(2)  p(l)  .  p(2)  ,(1)  .  v<=> 

zb  zb  ^  b  ^  b  .  zb  zb 


Ob 


Substituting  Eqs.  (3.119)  and  (3.120)  into  the  two  above  equations  yields  the 
following  expressions: 


id) ,  ,(2) 
zb  zb 
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and 


Al)  A2) 
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f, 


Therefore,  it  follows 
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b.  Solenoidal  components  of  the  beam  current  density. 
We  have 


It  is  known  that 


(Sb)*- 


l<3) 

zb 
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JF.hz 


Therefore,  we  have 


It  had  been  shown  before  [Eq.  (3-95)]  that 


W  2^0 


JF^^hz 


and 
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By  using  these  relations,  we  have 


On  the  other  hand,  by  using  Eqs.  (3. 105)  and  (3.117)^  the  following 
expression  is  obtained: 


v(3)  ^ 

I  I 
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v(i)  + 

I _ _I 

2V^ 


and  from  Eqs.  (3.110)  and  (3.1l6)  we  obtain 
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Substituting  the  latter  two  expressions  into  Eq.  (B.2),  we  obtain 
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Finedly,  the  total  electron  beam  current  density  is  given  by  the  sum  of 
the  right-hand  sides  of  Eqs.  (B.l)  and  (B.3)>  that  is, 


i  ,  (i  ,  )n.s.  (i  ,  )s. 
zb  _  ^  zb^  ^  zb^ 


✓ 

✓  \ 

_  _ 

\ 

aM 

( 

-j|z 

s  - < 

J  1  —  1  sin(hz)  + 

e  •  cos(Az)  -  cos(hz) 

2 

vJ 

- 

The  toted  volume  current,  at  the  output  grids,  (z  =  L)  ,  is  cal¬ 
culated  by  averaging  over  the  cross-section  of  the  beam  the  quantity 
(i^^i^)  and  by  mvdtiplying  the  result  by  the  beam  coupling  coefficient 
M  .  We  have 


+ 


cos(fJ^L)  -  cos(hL) 


Figure  U  shows 


as  a  function  of  (P^a)  . 

2.  Plasma  Electrons  (in  Region  l) 

The  nonsolenoidal  cooiponent  of  the  plasma  current  density  is 
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But,  it  is  known  that  the  lineeurized  longitudinal  plasma  current  densities 
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Then,  it  follows  that 
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By  substituting  Eq. .(3.112)  into  the  two  above  equations,  we  get 
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Therefore,  it  follows  that 
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For  solenoldal  component  of  the  plasma  current  density,  we  have 
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But  we  know  [Eqs.  (3*65)]  that 
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By  using  the  expressions  derived  in  the  text  for  these  potentials,  V'^ 
and  [Eqs.  (3.IIO)],  the  solenoidal  part  of  the  plasma  volume 

current  density  turns  out  to  be  given  by  the  following  expression: 
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The  total  plasma  volume  current  density,  at  emy  point  in  the  plasma-beam 
interaction  region,  distant  z  from  the  input  grids  (O  <  z  <  L)  ,  is 
therefore 
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APPENDIX  C 

CALCULATION  OF  THE  SURFACE  CURRENTS 

To  calculate  the  surface  currents  per  unit  of  circuntference,  at 
r  ss  a  ,  we  have  to  determine  the  surface  currents  originated  by  the 
electrons  of  the  beam  in  the  two  propagating  waves  modes.  We  have,  then, 
nonsolenoidal  surface  currents,  which  are  given  by 


h  dr  \  2V, 


r=a 


and 


(C.l) 


^  1-  dr  V2V„  , 

0  \  0  /rea 


.(1,2) 


Froa  the  condition  that  (Aa)  must  vanish  at  z  *  0  ,  we  obtain 
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It  follows,  from  Eqs.  (C.l)  and  (3.118),  that  the  right-hand  side  of  the 
above  equation  can  be  written  as 
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(i  ^)n.s.  aM  / A 

— rr 


80^  I  lo^e^) 


h  \  I,(P  a) 

-  J  I  —  j  sin(hz)  ® 


+  ... 


For  the  solenoidal  components,  we  have 

r(3)' 


sb 


d  /V' 


[cos(hz)  -  1] 
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(C.2) 
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and 


i(M 

sb 


At  z  =  0  ,  we 


r=a 


dr\2Vo 


rssa 


must  have  =  0  ,  from  which  we  obtain 
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Therefore,  from  Eqs.  (C,2)  and  (C.3),  the  total  surface  current  density, 
at  r  *  a  due  to  the  electrons  of  the  beam,  is 


sln(hz)  sln(Az) 

-  e  2 

®  h  ®  A 

—  w 

The  total  surface  current,  at  z  »  L  (output  grids),  is  determined  by 
averaging  the  surfewe  current  density  over  the  cross-section  of  the  beam 
and  multiplying  the  result  by  the  beam  coupling  coefficient  M  . 
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The  result  is 


sin(Az)  -j|z  sin(hz) 
-  e  -  - 
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APPENDIX  D 


CALCULATION  FOR  b/bt  (l^  ) 

In  this  appendix,  we  will  go  through  some  purely  mathematical 

manipulations  which  would  lead  to  an  expression  for  b/bt  (  1^  )  in 

D,  2 

terms  of  the  fundamental  and  second-harmonic  components  of  the  electric 
field.  From  Eq.  (^.35)^  we  have 
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It  can  be  shown,  however,  from  Eq  (4.15),  that 
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and  from  Eqs.  (4.13)  and  (4.15)>  the  following  expression  results: 
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From  the  same  eqxiations  (4.13)  and  (4.15),  we  obtain 
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It  follows,  therefore,  that 
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To  calculate  (a/at)(i^  we  make  use  of  Eq.  (4.3^)  whose  partial 

time  derivative  is  ^ 
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From  Eq.  (4.11)  we  derive  the  result 
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From  Eqs.  and  (4.l4)  we  obtain 
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from  which  we  derive  the  result 
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It  follows,  then,  that 
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APPENDIX  E 

DERIVATION  OF  THE  PARTICULAR  SOLUTION  OF  THE  DIFFERENTIAL  EQUATION  (4.43) 

The  inhomogeneous  term  of  the  differential  equation  can  be  written 
in  the  following  form: 
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Consider  the  part  in  the  above  expression,  which  corresponds  to  the  phase- 
2  iX 

factor  e  ,  where  X  =  a>t-P  z  .  We  have 
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The  complementary  function  is  written 

D2(z,t)  =  (k^ 

where  and  are  arbitraury  constants. 
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According  to  the  method  of  variation  of  parameters  to  determine  a 
particular  integral  for  the  differential  equation,  we  assume  that 
and  Bg  are  functions  of  z  ,  which  should  now  be  looked  upon  as  two 
new  dependent  variables. 

The  following  condition  is  imposed  on  A^Cz)  and  B2(z)  : 

dA^Cz)  Jh^z  ^^{z)  -jhgZ 

-  e  +  -  e  =  0  . 

dz  dz 

From  this  system  of  equations, 
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dz  dz 


Inhomogeneous  Term 


dAg/dz  and  dBg/dz  are  determined.  We  can  then  find  expressions  for 
Ag  and  Bg  by  direct  integration. 

The  results  eure 
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It  follows  thon  that  tho  particular  integral  corresponding  to  that  part 

2  iX 

of  the  inhomogeneouG  term  with  the  phase  factor  e  is 
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APPENDIX  F 

DERIVATION  OF  SECOND  HARMONIC  VELOCITIES 


The  second  harmonic  velocity  is  calculated  from 
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where  P  is  defined  as  being  given  by 
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From  Eqs.  (F.2)  and  (F.3),  the  following  expression  is  obtained  for  : 
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It  follows,  from  (F. 4),  that 
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By  using  the  expression  for  P  ,  we  have 
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By  substituting  this  above  expression  and  (F.5)  in  (F.l),  we  obtain 
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Ithaca,  New  York 

ATTN:  Prof.  G.  C.  Dalman 


University  of  Southern  California  1 

University  Park 

Los  Angeles  7,  California 

ATTN:  Z.  A.  Kaprielian 

Associate  Professor  of  Electrical  Engineering 

Dr.  A.  L.  Cullen  1 

Department  of  Electrical  Engineering 
University  of  Sheffield 
Sheffield  1,  England 
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No.  of  Copies 


Code 

u  358 

U  359 

u  391 

u  430 

G  70 

K  160 

I  13 

I  53 

I  96 

I  237 


Organization 

University  of  Arizona  1 

Tucson  25,  Arizona 

ATTN:  Prof.  Donald  C.  Stinson 

Department  of  Electrical  Engineering 

The  University  of  British  Colxanbia  1 

Department  of  Electrical  Engineering 
Vancouver  8,  B.  C.  Canada 
ATTN:  G.  B.  Walker 

Microwave  Laboratory 

University  of  Illinois  1 

Electrical  Engineering  Research  Laboratory 
Urbana,  Illinois 
ATTN:  Professor  P.  D.  Coleman 
Ultramicrowave  Group 

Chalmers  University  of  Technology  1 

Glbraltargaton  5  G 

Gothenburg,  Sweden 

ATTN:  Mr.  H.  Wilhelmsson 

Research  Laboratory  of  Electronics 

» 

Advisory  Group  on  Electron  Devices  (AGED)  2 

Office  of  the  Director  of  Defense  Res  and  Eng. 

346  Broadway,  8th  Floor 
New  York  13,  New  York 

U.  S.  Naval  Reseeu”ch  Laboratory  1 

Washington  25,  D.  C. 

ATTN:  Mr.  F.  J.  Liberatore,  Code  7420 

Bell  Telephone  Laboratories,  Inc.  2 

Whippany  Laboratory 

Whippsuiy,  New  Jersey 

ATTN:  Technical  Information  Library 

Hughes  Aircrsift  Coo^any  1 

Florence  Avenue  and  Teale  Street 
Culver  City,  Calif. 

ATTN:  Documents  Section 

Research  and  Development  Library 

Sandia  Corporation  1 

Sandia  Base,  P.  0.  Box  580O 

Albuquerque,  New  Mexico 

ATTN:  Mrs.  B.  R.  Allen,  Librarian 

Technical  Library  1 

G.  E.  TWT  Product  Section 

601  California  Avenue 

Palo  Alto,  California 

ATTN:  Verna  Van  Velzer,  Librarian 


-  5  - 


No.  of  Copies 


Code 
I  260 

I  266 

I  297 

I  305 

I  306 

I  309 

I  310 

I  312 

I  366 


Organization 

Sylvania  Elec.  Prod.  Inc.  1 

Electronic  Defense  Laboratory 

P.  0.  Box  205,  Mountain  View,  Calif. 

ATTN:  Library 

ITT  Federal  Laboratories  1 

Technical  Library 
500  V/ashington  Avenue 
Nutley  10,  New  Jersey 

Sperry  Gyroscope  Company  1 

Division  of  Sperry  Rand  Corporation 
Great  Neck,  N.  Y. 

ATTN:  Mrs.  Florence  W.  Turnbull 
Engineering  Librarian 

General  Electric  Company  1 

Power  Tube  Dept.  Electronic  Con^onents  Division 
Building  269,  Room  205 
One  River  Road 
Schenectady  5>  New  York 

General  Telephone  and  Electronics  Laboratories,  Inc.  1 
Bayside  Laboratories 
Bayslde  6,  New  York 
ATTN:  Mr.  D.  Lsaare 

Manager,  Project  Aim. 

Litton  Industries,  Inc.  1 

960  Industrial  Road 

San  Carlos,  California 

ATTN:  Document  Custodian,  Engr.  Dept. 

Veirlan  Associates  1 

611  Hansen  Way 

Palo  Alto,  California 

ATTN:  Dr.  Richard  B.  Nelson 

STL  Technical  Libreuy  1 

Document  Acquisitions 

Space  Technology  Laboratories,  Inc. 

P.  0.  Box  95001 

Los  Angeles  45,  California 

Radio  Corporation  of  America  1 

Defense  Electronic  Products 
Camden,  New  Jersey 

ATTN:  Mr.  S.  Schach,  Building  10-5 

Standards  Engineering,  Section  577 


No,  of  Copies 


Code 
I  367 

I  370 

I  382 

I  384 

I  435 

I  450 

I  547 

I  562 

I  577 

I  594 


Organization 

Stanford  Research  Institute  1 

Document  Center 

Menlo  Park,  California 

ATTN:  Acquisitions 

General  Telephone  and  Electronics  Labs.,  Inc,  1 

Bayside,  ..ew  York 

ATTN:  Dr.  T.  G,  PoLanyi,  Head  Thermioaics  Branch 

RCA  Laboratories  1 

David  Sarnoff  Research  Center 
Princeton,  New  Jersey 
ATTN:  Dr.  Harwick  Johnson 

Bell  Telephone  Laboratories  2 

Murray  Hill  Laboratory 
Murray  Hill,  New  Jersey 
ATTN:  Dr.  J.  R.  Pierce 

General  Electric  Company  1 

P.  0.  Box  1088 

Schenectady,  New  York 

ATTN:  Mr.  E.  D.  McArthur,  Manager 

Superpower  Microwave  Tube  Laboratory 

Bell  Telephone  Laboratories  1 

Murray  Hill,  New  Jersey 
ATTN:  J.  W.  Fitzwilliam,  Director 
Electron  Tube  Development 

The  Rand  Corporation  1 

1700  Main  Street 

Santa  Monica,  California 

ATTN:  Technical  Librarian 

Philips  Laboratories  1 

A  Division  of  North  American  Philips  Co., Inc. 
Irvington  on  Hudson,  New  York 
ATTN:  William  P.  Arnett 
Security  Officer 

Raytheon  Manufacturing  Con^iany  1 

520  Winter  Street 
Waltham,  Massachusetts 
ATTN:  Mr.  0.  T.  Fundingsland 

Research  Technology  Associates,  Inc.  1 

100  Ladge  Drive 

Electronic  Park  at  Avon 

Avon,  Massachusetts 

ATTN:  J.  Babakian 


-  7  - 


Code 
I  595 

I  666 

I  756 

I  759 

I  760 

I  761 

I  762 

I  763 

I  764 

I  765 

I  766 


Organization 

Giannlnni  Research 
Santa  Ana,  California 
ATTN:  J.  K.  Hagele 

Technical  Librarian 

Ramo-Wooldrldge 

A  Division  of  Thompson  Ramo-Wooldrldge  Inc. 

8433  Fallbrook  Avenue 

Canoga  Park,  California 

ATTN:  Technical  Information  Services 

Varian  Associates 
611  Hansen  Way 
Palo  Alto,  California 
ATTN:  Mr.  C.  W.  McClelland 

Technical  Publications  Manager 

Stanford  Research  Institute 
Menlo  Peurk,  California 
ATTN:  Mr.  C.  J.  Cook 

General  Atomic,  Div.  of  General  dynamics  Corp. 
P.  0.  Box  608,  San  Diego,  California 
ATTN:  Mr.  M.  Rosenbluth 

Linfield  Research  Institute 

McMinnville,  Oregon 

ATTN:  Dr.  V/.  P.  Dyke,  Director 

Coliimbia  Radiation  Laboratory 
538  West  120th  Street 
New  York  27,  New  York 

Sperry  Gyroscope  Conqpany 

Engineering  Library 

Mall  Station  C-39 

Great  Neck,  Long  Island,  New  York 

Watkins -Johns on  Conpany 
3333  Hillvlew  Avenue 
Palo  Alto,  California 
ATTN:  Dr.  H.  R.  Johnson 

Westinghouse  Electric  Corporation 
Friendship  International  Airport 
Box  746.,  Baltimore  3,  Maryland 
ATTN:  G.  Ross  Kilgore,  Manager 

y^plied  Research  Depeurtment 
Baltimore  Laboratory 

Services  Electronic  Research  Laboratories 
Baldock,  Herts,  England 
ATTN:  Dr.  Boot 


No.  of  Copies 
1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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Code 
I  767 

I  768 

I  769 

I  770 

I  999 

I  771 

I  772 

I  9^^ 

U  2 

U  21 

U  22 


Organization 

Standard  Telephone  Laboratories 
Harlow,  Essex,  England 
ATTN:  Dr.  E.  A.  Ash 

Elt  e  1-Mc  Cullo\agh ,  Inc . 

798  San  Mateo  Avenue 
San  Bruno,  California 
ATTN:  Librarian 

Hewlett-Packard  Company 
275  Page  Mill  Road 
Palo  Alto,  California 

Hughes  Aircraft  Company 

Research  and  Development  Laboratories 

Culver  City,  California 

ATTN:  L,  M.  Field 

Raytheon  Company 

Box  171,  Way  land.  Mass. 

ATTN:  E.  Rolfe 

Advanced  Kinetics  Inc. 

P.  0.  Box  1803 

Newport  Beach,  California 

ATTN:  Dr.  R.  Wanlek 

RCA  Laboratories 
Princeton,  New  Jersey 
ATTN:  Harwell  Johnson 

Bomac  Laboratories,  Inc. 

8  Salem  Road 
Beverly,  Massachusetts 
ATTN:  Dr.  Arthur  McCoubrey,  Manager 
Research  and  Development 

California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
Pasadena  4,  California 
ATTN:  Documents  Library 

The  Johns  Hopkins  University 
Radiation  Laboratory 
1315  St.  Paul  Street 
Baltimore  2,  Maryland 
ATTN:  Librarian 

The  Johns  Hopkins  University 
Department  of  Physics 
Homewood  Campus 
Baltimore  I8,  Maryland 
ATTN:  Dr.  Donald  E.  Kerr 


No.  of  Copies 
1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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No.  of  Copies 


Code 
U  26 

U  4o 

U  42 

U  59 


U  79 


U  100 

U  102 


U  107 


Organization 

Massachusetts  Institute  of  Technology  1 

Lincoln  Laboratory 
P.  0.  Box  73 

Lexington  73,  Massachusetts 
ATTN:  Mary  A,  Granese,  Librarian 

New  York  University  1 

Department  of  Physics 

College  of  Arts  and  Sciences 

Washington  Square,  New  York  3,  New  York 

ATTN:  Professor  J.  H.  Rohrbaugh 

The  Ohio  State  University  1 

2024  Nell  Avenue 

Colianbus  10,  Ohio 

ATTN:  Prof.  E.  M.  Boone 

Department  of  Electrical  Engineering 


Library  1 

Georgia  Technology  Research  Institute, 

Engineering  Experiment  Station 
722  Cherry  Street,  N.  W. 

Atlanta,  Georgia 

ATTN:  Mrs.  J.  H.  Crosland,  Librarian 

University  of  Michigan  1 

Engineering  Research  Institute 

Radiation  Laboratory 

AITN:  Prof.  K.  M.  Siegel 

912  N.  Main  Street 

Ann  Arbor,  Michigan 

University  of  Ceilifornia  1 

Electronics  Research  Lab. 

332  Cory  Hall,  Berkeley  4,  Calif. 

ATTN:  J.  R.  Whinnery 

Harvard  University  1 

Technical  Reports  Collection 

Gordon  McKay  Library 

303A  Pierce  Hall 

Oxford  Street 

Cambridge  38,  Massachusetts 
ATTN:  Librarian 

University  of  Michigan  1 

Electronic  Defense  Grovq) 

Engineering  Research  Institute 

Ann  Arbor,  Michigan 

ATTN:  J.  A.  Boyd,  Supervisor 


-  10  - 


Code 
U  Il4 

u  150 

u  168 

u  169 

u  209 

u  228 

u  237 

U  240 


U  284 


Organization 

University  of  California 
390  Cory  Hall 
Berkeley  4,  Calif. 

ATTN:  Dr.  Charles  Susskind 

University  of  Washington 
Department  of  Electrical  Engineering 
Seattle  5^  Washington 
ATTN:  Mr.  A.  E.  Harrison 

The  George  Washington  University 
Department  of  Electrical  Engineering 
School  of  Engineering 
Washington  6,  D.  C. 

ATTN:  Dr.  N.  T,  Grisamore,  Exec.  Officer 


No.  of  Copies 
1 


1 


1 


1 


Illinois  Institute  of  Technology  1 

3301  S.  Dearborn  Street 
Chicago  16,  Illinois 
ATTN:  Dr.  George  I.  Cohn 


New  York  University  1 

College  of  Engineering 

346  Broadway 

New  York,  New  York 

ATTN:  Mr.  L.  S.  Schwartz,  Research  Division 


University  of  Kansas  1 

ElectriceJ.  Engineering  Depeurtment 
Lawre  nee ,  Kans  as 
ATTN:  Dr.  H.  Unz 

Polytechnic  Institute  of  Brooklyn  1 

Microwave  Research  Institute 

55  Johnson  Street 

Brooklyn,  New  York 

ATTN:  Dr.  N.  Marcuvitz 

Illinois  Institute  of  Technology  1 

Technology  Center 

Department  of  Electrical  Engineering 
Chicago  16,  Illinois 
ATTN:  Paul  C.  Yuen 

Electronics  Research  Laboratory 


Dr.  J.  T.  Senise  1 

Instituto  Tecnologico  de  Aeronaut ica 
Sao  Jose  dos  Caoqpos 
Sao  Paulo,  Brasil 


11  - 


Code 
U  288 

u  292 

u  294 

u  308 

u  320 

u  321 

u  322 

u  323 

U  324 

u  325 

u  326 


Organization 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
55  Johnson  Street 
Brooklyn  1,  New  York 

University  of  Maryland 
College  Park,  Maryland 
ATTN:  Dr.  J.  M.  Burgers 

University  of  Illinois 

Electrical  Engineering  Reseeurch  Laboratory 

Urbana,  Illinois 

ATTN:  Dr.  A.  A.  Dougal 

Brandels  University 
Waltham,  Massachusetts 
ATTN:  Dr.  E.  P.  Gross 

California  Institute  of  Technology 
1201  East  California  Street 
Pasadena,  California 
ATTN:  Prof.  R.  W.  Gould 

University  of  California 
Radiation  Laboratory 
Livermore,  California 
ATTN:  N.  Chrlstofllos 

University  of  California 
Radiation  Laboratory 
Livermore,  California 
ATTN:  S.  A.  Colgate 

University  of  California 
Radiation  Laboratory 
Livermore,  California 
ATTN:  R.  F.  Post 

Princeton  University 
Princeton,  New  Jersey 
AITN:  L.  Spltzer 

Project  Matterhorn 

Massachusetts  Institute  of  Technology 
77  Massachusetts  Avenue 
Cambridge,  Massachusetts 
ATTN:  W.  P.  Allis 

University  of  California 
Electrical  Engineering  Department 
Berkeley  4,  California 
ATTN:  Prof.  J.  R.  Singer 


No.  of  Copies 
1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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No.  of  Copies 


Code 

u  327 

u  328 


u  329 

u  330 

u  331 

u  332 

U  333 

u  334 

U  335 

U  337 


Organization 

University  of  California  1 

Radiation  Laboratoiy 
Berkeley,  California 
ATTN:  Dr,  R,  K.  Wakerling 

Information  Division,  Building  ^0,  Room  128 

University  of  Chicago  1 

Institute  of  Air  Weapons  Research 

Museum  of  Science  and  Industry 

Chicago  37,  Illinois 

ATTN:  Mrs.  Norma  Miller, 

Technical  Librarian 

University  of  Florida  1 

Department  of  Electrical  Engineering 

Gainesville,  Florida 
ATTN:  Prof,  W,  E.  Lear 

University  of  Illinois  1 

Control  Systems  Laboratories 

Urbana,  Illinois 

ATTN:  Prof.  Daniel  Alpert 

University  of  Illinois  1 

Department  of  Electrical  Engineering 

Urbana,  Illinois 
ATTN:  Prof.  L.  Goldstein 

University  of  Illinois  1 

Department  of  Physics 

Urbana,  Illinois 

ATTN:  Dr.  John  Bardeen 

Johns  Hopkins  University  1 

Radiation  Laboratory 

1315  St.  Paul  Street 

Baltimore  2,  Maryland 

ATTN:  J.  M.  Minkowski 

Iowa  State  University  1 

Physics  Department 
Iowa  City,  Iowa 

ATTN:  Professor  Frank  MacDonald 

University  of  Colorado  1 

Department  of  Electrical  Engineering 
Boulder,  Colorado 
ATTN:  Professor  W.  G.  Worcester 

McMurray  College  1 

Department  of  Physics 

Abilene,  Texas 

ATTN:  Dr.  Virgil  E.  Bottom 


-  13  - 


Code 

Organization 

No.  of  Copies 

U  330 

University  of  Michigan 

3506  East  Engineering  Building 

Ann  Arbor,  Michigan 

ATTN:  Electron  Tube  Laboratory 

1 

U  339 

University  of  Minnesota 

Institute  of  Technology 

Department  of  Electrical  Engineering 
Minneapolis,  Minnesota 

ATTN:  Prof.  A.  Van  der  Ziel 

1 

u  3^0 

Ohio  University 

College  of  Applied  Science 

Athens,  Ohio 

ATTN:  D.  B.  Green 

1 

u  3^1 

Oregon  State  College 

Department  of  Electrical  Engineering 
Corvallis,  Oregon 

ATTN:  H.  J.  Oorthi^ys 

1 

u  3U2 

Princeton  University 

Department  of  Electrical  Engineering 
Princeton,  New  Jersey 

1 

u  3^3 

Purdue  University 

Research  Library 

Lafayette,  Indiana 

ATTN:  Electrical  Engineering  Department 

1 

u  344 

Rensselaer  Polytechnic  Institute 

Office  of  the  Llbrarisui 

Troy,  New  York 

1 

u  345 

Scientific  Attache 

Swedish  Elnbassy 

2249  R  Street,  N.  W. 

WMhington  8,  D.  C, 

1 

U  346 

Rutgers  University 

Physics  Department 

Newark  2,  New  Jersey 

ATTN:  Dr.  Charles  Pine 

1 

u  347 

University  of  Texas 

Military  Physics  Research  Laboratory 

Box  8036,  University  Station 

Austin,  Texas 

1 

U  348 

University  of  Utah 

Electrical  Engineering  Department 

Salt  Lake  City,  Utah 

ATTN:  Richard  W.  Grow 

1 

-  l4  - 

No«  of  Copies 


Code 
U  350 

U  351 

U  352 

u  353 

U  35*^ 

U  355 

U  356 

U  436 


Organization 


University  of  Puerto  Rico  1 

College  of  Agriculture  and  Mechanical  Arts 
Mayaguez,  Puerto  Rico 
ATTN:  Dr.  Braulio  Dueno 

The  Royal  Institute  of  Technology  1 

Stockholm  JO,  Sweden 
ATTN:  Dr.  B.  Agdur 

Drexel  Institute  of  Technology  1 

Department  of  Electrical  Engineering 
Philadelphia  4,  Pennsylvania 
ATTN:  F.  B.  Haynes 

Forsvarets  Forskningsinstitutt  1 

Avdeling  for  Radar 
Bergen,  Norway 

ATTN:  WADD  (WCOSR,  Mr.  Knutson) 

Wrlght-Patterson  AFB,  Ohio 

Massachusetts  Institute  of  Technology  1 

Research  Laboratory  of  Electronics 
Cambridge  39,  Mass. 

ATTN:  L.  D.  Smullin 

Massachusetts  Institute  of  Technology  1 

Research  Laboratory  of  Electronics 
Cambridge  39,  Mass, 

ATTN:  S.  C,  Brown 

University  of  Chicago  1 

Midway  Laboratories 

6220  S.  Drexel  Avenue 

Chicago,  Illinois 

ATTN:  P.  J.  Dickerman 

University  of  Mississippi  1 

University,  Mississippi 
ATTN:  I'lr.  Thomas  Tullos 

Hq.  AFCRL,  Office  of  Aerospace  Research  (CRRC)  4 

L.  G.  Hanscom  Field, 

Bedford,  Mass, 


-  15  - 


